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ABSTRACT 


A  detailed  analysis  is  presented  for  a  circular  array  of  Beverage 
antennas  over  a  conducting  earth  Equations  are  derived  for  azimuth  and 
elevation  patterns  for  any  number  of  similar  antennas  in  the  array  taken 
any  number  at  a  time,  spaced  in  any  manner,  and  summed  in  any  arbi¬ 
trary  manner  with  or  without  phasing,  ^f  particular  interest  for  HF 
direction  fijidii^  applications  are  sector^rrays  of  10  to  20  antennas 
located  and  operated  between  l^nd  40  me  Simple  summing 

•{nonphasodf  of  typical  sector  arrays  i^oduces  patterns  which  have 
greatly  reduced  beamwidth  compar^  to  single  element  patterns  At 
10  me,  a  25  nrieter  long  single  e^i^ment,  one  meter  high,  has  a  Kdb 
beamwidth  of  about  78®,  whU^e^  circular  sector  array  of  21,^^fmilar 
elements  2®  apart  has  a^jx<amwidth  of  about  IS^witliout^asing  Ele¬ 
vation  patterns  element,  1  a rj^ .sec-tor  arrays  show  similarly 

narrow  azinmth’^tterns  .up  to  "eTevation  angles  of  70®  at  10  me  for  an 
array  dj^rt^ietej:  oFTess  than  900  feet. 

-i  he  analysis  provides  both  ground  wave  and  sky  wave  patterns 
For  sky  wave  patterns  the  polarization  of  the  incident  field  may  be 
linear  or  elliptical  iany-conditioii  between  vertical  aiTd  liui'lzunial  ptTl'dri- 
-aation)-.  Calculations  show  the  sector  arrays  considered  .so  far  to  have 
low  polarization  error.  Polarization  error  is  further  reduced  as  beam- 
width  is  reduced.  _ 


The  analysiXalso  provides  antenna  impedance,  antenna  line 
constants,  effective  \eight,  wave  tilt  angle  and  various  other  param¬ 
eters.  Calculated  pe^ormance  shows  theantennas  to  have  impedances 
and  patterns  normally Vssociated  with  frequency  independent  antennas. | 
The  results  show  that  D/F  performance  over  a  100  to  1  frequency  ranc^ 
extending  as  low  as  1  5  me,  should  be  obtained  in  a  single  array  incluaing 
good  azimuth  patterns  at  high  elevation  angles.  Furthermore,  the  sirn- 
plicity  of  the  antenna  element  design  permits  what  is  probably  Jhe^ -lowest 
cost  wide  aperture  direction  findex  antenna  array  yel^esigh^  for  the 
MF,  HF,  and^H^  fre-qqency  r^^e^T^  ^  ^ 

"S^^^alculated  '^■nH  effective 

heights  indicate  that  the  summed  sector  arrays  provide  an  adequate  sen¬ 
sitivity  with  resulting  effective  heights  ranging  from  a  few  meters  to  the 
vicinity  of  100  meter  s  .Bearing  sensitivity  is  probably  as  good  as  present 
day  instrumentation  cajnTirfitlze  either  in  a  scanning  or  fixed  beam  mode. 
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It  is  also  apparent  from  the  analysis  that  the  performance  of  the 
antenna  array  is  not  highly  dependent  on  a  specific  earth  conductivity  as 
in  the  case  of  some  D/F  antennas.  Reception  is  improved  by  an  increased 
wave  tilt  angle  usually  associated  with  poorly  conducting  earth.  Thus  the 
system  is  probably  more  adaptable  to  varying  site  terrain  than  is  the  usual 
D/F  system. 

This  is  Part  IV  of  a  four-part  report  and  is  concerned  primarily 
with  theory.  Part  II  is  concerned  with  system  design  and  performance 
for  long  range  HF  D/F  while  Part  III  is  a  detailed  report  on  measured 
accuracy  of  a  900-ft  diameter  array  of  180  antennas.  Part  I  was  submitted 
in  1962.  All  theoretical  background  material  necessary  to  the  analysis  has 
been  reproduced  in  the  Appendix 
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It  was  required  that  a  procedure  be  developed  for  deriving 
optimum  circular  arrays  of  Beverage  antennas  for  high  frequency  direction 
finding.  Criteria  for  an  optimum  array  are  bearing  accuracy,  a  narrow 
beam  azimuth  pattern  and  a  favorable  elevation  pattern  for  sky  waves. 

Although  bearing  accuracies  will  not  be  calculated  directly  in  the  analysis 
which  follows,  it  will  be  assumed  that  if  the  azimuth  pattern,  formed  by  a 
wide  aperture  distribution  of  antenna  elements,  cam  be  made  more  directive, 
then  the  resulting  bearing  information  will  move  in  the  direction  of  improved 
accuracy.  The  conditions  for  favorable  elevation  patterns  and  maximum 
effective  height  will  also  be  determined.  The  analysis  will  also  permit 
estimates  of  polarization  error. 

The  principal  variables  to  be  optimized  are  the  number  of 
antennas,  the  manner  of  summing,  antenna  spacing  geometry,  antenna 
location  or  overall  array  size,  antenna  length  and  height  above  ground. 

Other  variables  such  as  the  termination  geometry  and  eintenna  conductor 
geometry  will  usually  be  determined  more  on  the  basis  of  practical 
requirements  rather  than  on  optimization  of  electrical  characteristics. 

Certain  other  variables  will  not  be  considered  in  this  analysis,  such  as 
the  mutual  coupling  between  antennas,  it  being  assumed  that  acceptable 
conditions  can  be  specified  from  experience.  (The  reference  of  Appendix  I 
is  sufficiently  general  to  treat  mutual  impedance  theoretically  if  this  should 
be  of  interest  at  any  future  time.  ) 

The  procedure  followed  in  the  derivations  given  below  was  to 
assume  the  geometry  shown  in  Figures  1,  Z  and  3  and  to  assume  independent 
ground  and  sky  wave  components  arriving  at  the  same  azimuth.  The  sky 
wave  component  includes  both  vertical  and  horizontal  polarization  components 
which  may  be  specified  independently. 

The  total  field  strength  for  direct  and  reflected  rays  is  calculated 
at  a  point  on  the  antenna  wire,  with  phase  referred  to  the  origin.  The  exact 
reflection  coefficients  of  the  earth  for  both  polarizations  are  retained  in  the 
calculation  for  the  sky  waves,  and  the  earth  constants  as  they  determine  the 
wave  tilt  for  the  ground  wave. 

The  analysis  is  applied  to  the  kth  antenna  of  a  circular  array 
as  shown  in  Figure  1,  all  antennas  being  arranged  in  symmetrical  pairs 
about  the  x  axis.  All  antennas  are  assumed  to  have  similar  geometry  and 
electrical  characteristics.  The  earth  under  the  antennas  is  assumed  to  he 
constant  in  the  x  and  y  directions. 
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The  receiving  ond  current  of  the  kth  antenna  is  found  by 
integration  along  the  conductor.  The  currents  for  all  antennas  are  then 
summed  and  assigned  arbitrary  amplitude  and  phase  constants  such  as 
would  correspond  to  insertion  of  a  goniometer  or  special  phasing  networks 
for  fixed  beams. 

Analytical  methods  of  summing  antenna  element  currents  to 
determine  the  complete  array  patterns  are  various,  such  as  the  use  of 
Fourier  transforms.  Most  of  these  methods  are  useful  on  linear  arrays 
particularly  with  omnidirectional  elements.  In  the  case  of  radially  oriented 
end  fire  elements  in  a  circular  array,  no  significantly  useful  short  cut 
methods  have  been  found  which  allow  clear*cut  trading  between  beamwidth 
and  side  lobe  level  (for  instance).  It  has,  therefore,  been  decided  to  use  a 
digital  computer  to  sum  the  antenna  currents  derived.  For  the  present, 

"trial  and  error"  pattern  calculations  are  being  used;  that  is, most  calculations 
are  based  on  assumed  antenna  arrays,  with  substitution  of  measured 
electrical  characteristics  which  are  appropriate  to  the  local  site.  The 
resulting  patterns  are  then  evaluated  for  quality.  This  method  has  been 
found  to  be  quite  satisfactory. 

A  large  number  of  calculated  patterns  are  presented  in  this 
report  (Part  IV)  .  *At  some  later  date  it  may  be  advantageous  to  use  more 
advanced  programming  methods  to  design  arrays  for  specific  optimum 
conditions,  i.e.,  specified  side  lobe  level,  bandwidth,  etc.  At  the  present 
time,  however,  the  problem  to  be  solved  does  not  require  this  degree  of 
sophistication.  Some  of  the  more  interesting  patterns  calculated  have  been 
experimentally  checked  by  field  measurements.  A  few  examples  of  these 
checks  have  been  included  in  Part  II  of  this  report. 

2.  Symbol  Definitions  and  Some  Assumptions 

The  following  symbols  and  definitions  will  be  adopted  in 
addition  to  those  of  Figures  1.  2,  3  and  9.  All  parameters  are  in  MKS 
units  except  where  noted.  (The  use  of  emu  units  has  been  retained  in 
certain  parts  of  the  calculation  which  are  based  on  early  references). 


’i' Although  detailed  analysis  and  interpretation  of  the  patterns  has  been  placed 
in  Part  II. 


FIGURE  I 

CIRCULAR  ARRAY  OF  BEVERAGE  ANTENNAS.* 
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FIGURE  Z 

DETAIL  OF  THE  GROUND  WAVE  NEAR  THE  ANTENNA. 
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FIGURE  3 

DETAIL  OF  THE  EARTH  REFLECTION  NEAR  THE  ANTENNA. 
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TABLE  1.  SYMBOL  DEFINITIONS 


P 

Po 

y 

6 


c 

Cw 

n 

e 

\ 

FA 

^^o 

e 


PM-) 


rg 


vs 


(t> 

(j)g 


A 

32 

B 

b 


=  Attenuation  constant  (nepers  per  meter) 

=  Phase  constant  =  2k/\  (radians  per  meter) 

=  Free  space  phase  constant 
=  Parameter  designating  an  arbitrary  antenna 
=  Propagation  constant  =  a  +  jP 
=  Wave  tilt  angle 
=  Permittivity 

=  Relative  dielectric  constant  of  the  earth 
=  Permittivity  of  free  space 
=  Angle  defined  in  Figure  3 
=  Amplitude  of  Wise's  correction  factor 
=  Phase  of  Wise's  correction  factor 
=  Plus  or  minus  azimuth  of  kth  antenna 
=  Angle  of  incidence 

=  Wavelength,  or  wavelength  on  the  antenna 
=  Free  space  wavelength 

=  Permeability,  also  an  integration  parameter  in 
equation  (44) 

=  Permeability  of  the  earth 
=  Permeability  of  free  space 

=  Phase  angle  representing  path  length  difference  between 
direct  ray  and  reflected  ray  at  Point  Pk 
=  Arbitrary  phase  constant  for  k(-)th  antenna 
=  Arbitrary  phase  constant  for  k(+)th  antenna 
=  Conductivity 
=  Conductivity  of  the  earth 

=  Arbitrary  phase  constant  for  horizontally  polarized 
sky  wave 

=  Arbitrary  phase  constant  for  ground  wave 
=  Arbitrary  phase  constant  for  vertically  polarized 
sky  wave 

=  -  Azimuth  of  signal  arrival 
=  -  Azimuth  of  arrival  of  Eg 
=  -  Azimuth  of  arrival  of  Eg 

=  Angle  of  elevation  above  xy  plane  of  direction  of 
arrival  of  Es 

=  ZlTf 

=  Parameter  in  Wise's  analysis 
=  Arbitrary  constant 
=  Arbitrary  constant 
=  Parameter  in  Wise's  analysis 
=  Radius  of  the  antenna  conductor 
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TABLE  1.  SYMBOL  DEFINITIONS  (Cont'd) 

C  =  Capacitance  per  unit  length 
c  =  Velocity  of  light  =  3  X  10®  meters/second 
Egp  =  Component  of  ground  wave  electric  field  parallel 
to  kth  antenna  at  Pk 

Ejjg  =  Horizontally  polarized  part  of  ground  wave  electric 
field  (perpendicular  to  direction  of  travel) 

Ehgp  =  Component  of  horizontally  polarized  ground  wave  electric 
field  parallel  to  kth  antenna  at  Pj^ 

Ejjg  =  Horizontally  polarized  part  of  sky  wave  electric  field 

(perpendicular  to  plane  of  incidence  and  direction  of  travel) 
Ep  =  Total  electric  field  parallel  to  kth  antenna  at  Pk 
Ej.g  =  Horizontally  polarized  part  of  ground  wave  which  is  also 

in  direction  of  propagation,  sometimes  called  the  wave  tilt 
component  or  the  radial  component  (lies  parallel  to  the 
xy  plane) 

Ej.gp  =  Component  of  radial  ground  wave  electric  field  parallel 
to  kth  antenna  at  Pk 

Egp  =  Component  of  sky  wave  electric  field  parallel  to  the  kth 
antenna  at  Pk 

Eyg  =  Vertically  polarized  part  of  ground  wave 
E^g  =  Vertically  polarized  part  of  sky  wave  (in  plane  of  incidence 
and  perpendicular  to  direction  of  travel) 
f  =  F  requency 

G  =  Conductance  per  unit  length 

Gk(-)  =  Arbitrary  amplitude  constant  for  k(-)th  antenna 
Gk(+)  =  Arbitrary  amplitude  constant  for  k(+)th  antenna 
h  =  Height  of  antenna  above  ground 

he  =  Effective  height 

h'  =  Parameter  in  equation  (44) 

I  =  Current 

Ik  =  Current  of  kth  pair  modified  by  summing  networks 

Iko  =  Current  from  the  kth  pair 

Ik(-)  ~  Current  on  k(-)th  antenna  at  Pk 

lk(-)o  “  Current  from  k(-)th  antenna 

lk(+)  ~  Current  on  k(+)th  antenna  at  Pk 

lk(+)o  ‘  Current  from  k(+)th  antenna 

j  =  Imaginary  unit 

K  =  A  constant  proportional  to  square  root  of  resistivity  as  in 
equation  ( 67) 

k(-)  =  Denoting  the  k(-)th  antenna 

k(+)  =  Denoting  the  k(+)th  antenna 

L  =  Length  of  the  Beverage  antenna  in  meters,  also  inductance 
per  unit  length  in  henries  per  meter 
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TABLE  1.  SYMBOL  DEFINITIONS  (Cont'd) 

t  -  First  direction  cosine 
ig  =  First  direction  cosine  for  Ng 
=  First  direction  cosine  for  Ng 
=  First  direction  cosine  for  rj 
il  =  First  direction  cosine  for  ■?£ 
m  =  Second  direction  cosine,  also  a  constant  integer 
mg  =  Second  direction  cosine  for  Ng 
mg  =  Second  direction  cosine  for  Ng 
mi  =  Second  direction  cosine  for  ?! 

m2  =  Second  direction  cosine  for  r2 

Ng  =  Unit  vector  in  direction  of  propagation  of  ground  wave 

Ng  =  Unit  vector  in  direction  of  propagation  of  sky  wave 

N  =  Summing  index  in  computer  program 
n  =  Velocity  ratio  or  subscript  denoting  the  nth  antenna, 
also  the  third  direction  cosine 
ng  =  Third  direction  cosine  for  Ng 

ng  =  Third  direction  cosine  for  Ng 

ni  =  Third  direction  cosine  for  ?! 

n2  =  Third  direction  cosine  for  r2 

P  =  Parameter  used  in  Carson's  analysis  as  in  equation  (43), 
emu  units 

P'  =  P  X  10-7  (P'  is  in  MKS  units) 

Q  =  Parameter  used  in  Carson's  analysis  as  in  equation  (43), 
emu  units 

Q'  =  Q  X  10“^  (Q'  is  in  MKS  units) 

R  =  Resistance  per  unit  length 

Rq  =  Real  part  of  Zq 

Rr  =  Radiation  resistance 

Rh  =  Complex  reflection  coefficient  of  the  earth  for  horizontal 
polarization 

Ry  =  Complex  reflection  coefficient  of  the  earth  for  vertical 
polarization 

r  =  Distance  from  origin  to  P^ 

rc  =  Carson's  parameter  defined  by  equation  (49)  in  emu 
Tew  “  Wise's  modification  of  Carson's  parameter  rc  in  emu 
rj^  =  Distance  to  Pj^  from  point  (Z  =  h,  y  =  0,  x  =  0) 

Tq  =  Inner  radius  of  the  array  circle 

r  =  Outer  radius  of  the  array  circle 

Tg  =  Distributed  loss  resistance  in  Carson's  analysis 

T  =  Vector  in  radial  direction  from  origin  to  P^ 

ri  =  Arbitrary  vector 

r2  =  Arbitrary  vector 
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TABLE  1.  SYMBOL  DEFINITIONS  (Cont'd) 

8  =  Wise's  correction  factor 

t  =  Time 

V  =  Voltage 

Vp  =  Voltage  induced  by  Ep 

V  =  Signal  velocity  on  the  antenna 

Xq  =  Imaginary  part  of  Zq 

Zq  =  Characteristic  impedance 

Zl,  =  Load  impedance,  terminating  impedance  of  the  line 
Zin  =  Impedance  locking  into  the  line  at  the  point  Pk 
Zj.  =  Sending  end  (receiver  end)  impedance  of  the  line 

Certain  other  symbols  and  equivalents  to  the  above  symbols 
are  given  in  Section  D  oj  this  report  where  the  digital  computer  program 
is  described. 


A  previous  analysis  by  A.  D.  Bailey^  considered  the  problem 
of  a  single  Beverage  antenna.  The  analysis  which  follows  is  similar  and 
extends  results  to  the  case  of  the  array  described  in  Figure  1. 

We  first  wish  to  establish  the  total  electric  field  intensity  at 
an  arbitrary  point  on  the  kth  antenna,  the  induced  current  will  then  be 
summed  over  the  length  of  the  antenna.  Proceeding  with  the  ground  wave 
Eg  arriving  at  the  azimuth  direction  -c|)g,  for  the  components  parallel  to  the 


antenna  at  Pk  we  have 

Ej-gp  =  Ej-gcos  (4)g  -  T]k) 

(1) 

Ehgp  ■  0 

(2) 

In  extreme  circumstances  such  as  a  very  close  target  trans¬ 
mitter,  or  very  unusual  values  for  the  earth  constants  CTg  and  Cg,  it  might 
occur  that  the  horizontally  polarized  ground  wave  E^g  is  not  equal  to  zero, 
hence  equation  (2)  is  an  approximation.  Nevertheless,  it  is  valid  under 
alniost  all  circumstances  of  practical  interest,  especially  when  a  vertically 
polarized  far  zone  transmitter  is  used.  The  case  of  the  antenna  considered 
as  a  source  yields  a  solution  which  bears  this  out  as  the  E^g  component 

1.  Bailey,  A.  D.  ,  Dyson,  J.D.,  and  Hayden,  E.  C.,  "S  udies  and  Investi¬ 
gations  Leading  to  the  Design  of  a  Radio  Direction  Finder  System  for  the 
MF-HF-VHF  Range,  '  '  Fourth  Quarterly  Report,  30  June  1962,  Contract 
No.  DA  36-093-SC- 87264,  EERL,  University  of  Illinois,  Urbana,  Illinois. 
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tangential  to  the  ground  "varies  as  the  inverse  square  of  the  distance  and 
is  of  no  significance  at  radio  frequencies^.  " 

Proceeding  with  the  sky  wave  there  are  four  components, 
two  direct  and  two  reflected  (or  two  vertical  and  two  horizontal),  all  arriving 
from  -  azimuth  (})g 

Egp  =  Ev8  cos  (4)s  -  nk)  sin  ijig  +  Ehs  sin  (<{>s  •  nk) 

+  Evs  Rve^^  cos  ((j>s  -  nk)  sin  +  Ehs  Rhe’^  sin  (4)3  -  nk) 

(3) 

There  is  no  radial  component  of  the  sky  wave  since  the  medium 
above  the  earth  is  nonconducting.  The  total  parallel  electric  field  at  the 
point  Pk  is  therefore 

Ep  =  ^gp  ^  ^sp  (^) 

where  Egp  is  the  sum  of  equations  (1)  and  (2). 

A  plane  earth  is  assumed  with  electromagnetic  constants  which 
range  within  reasonable  values  such  as 

Kg  "  Ko 
1  <  e  <  100 

o 

10"^  4  (Tg  4  5  mks 

It  has  also  been  assumed  that  the  transmission  efficiency  of  the 
line  is  reasonable,  that  is,  the  attenuation  constant  of  the  line  is  not  excessively 
large,  and  the  phase  velocity  on  the  line  is  somewhere  in  the  vicinity  of  the 
velocity  of  light  (at  least  .  3c). 

3.  Derivation  of  the  Total  Electric  Field  at  a  Point  on  the 

k t h  Antenn ^ 

From  inspection  of  Figure  1,  which  is  based  on  a  conventional 
spherical  coordinate  system,  the  direction  cosines  of  any  radial  direction 
from  the  origin,  defined  by  the  vector  r,  are  by  definition 


2.  Wait,  James  R.,  "Radiation  from  a  Ground  Antenna,  "  Canadian  Journal 
of  Technology,  Vol.  32,  No.  1.  January  1954.  pp.  1-9. 
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i  =  sin  U  sin  <> 

m  =  sin  d  cos  <}>  (5) 

n  =  cos  0 

The  cosine  of  the  angle  between  any  two  directions  defined  by 
unit  vectors  r  j  and  r  ,  is  therefore 

?!  ri  =  f  1  <2  +  ni  im2  njn^  (6) 

Now  consider  the  ground  wave  component  at  the  point  on 
the  k(-)th  antenna  as  indicated  in  Figures  1  and  2.  Note  that  the  planes  of 
constant  phase  are  tilted  forward  by  an  angle  6.  The  coordinates  of  the 
point  Pk  are 


X  =  I  r  I  cos  ^  cos  nk 
y  =  I  r  I  cos  ^  sin 
z  =  I r  I  sin  ^ 


(7) 


The  direction  cosines  of  the  unit  vector  Ng  normal  to  the  plane 
of  constant  phase  (wavefront)  of  the  ground  wave^  are 


fg  =  -cos  6  cos  (tJg 


mg  =  -cos  6  sin  4) 


g 


(8) 


Note  that  the  ground  wave  is  assumed  to  arrive  from  a  direction 
above  the  x-y  plane  by  the  amount  of  the  wave  tilt  angle  6.  The  details  of 
the  wave  tilt  phenomena  are  covered  briefly  in  a  later  section  and  are 

3.  In  general  Ej-g  and  Eyg  are  not  in  phase,  in  which  case  what  is  meant  is 
the  direction  of  the  minor  axis  of  the  ellipse  formed. 
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thoroughly  reported  in  the  literature^"^.  The  wave  tilt  angle  normally 
ranges  from  zero  (or  near  zero)  over  highly  conducting  surface's  to  some¬ 
thing  less  than  25*  over  most  dry  soil  surfaces.  Wait^  has  shown  that  this 
may  be  increased  over  a  stratified  ground.  Typical  values  are  given  in  a 
later  section  of  this  report. 


The  equation  of  the  plane  of  constant  phase  which  passes 
through  the  point  is  therefore  Ng  •  T  =  a  constant®.  Therefore 

Ng  •  r  =  f  g.x  +  mgY  +  OgZ  =  -  |r|(cos  ^  cos  6  cos  cos  <}>g 

+  cos  ^  cos  6  sin  sin  4)g  -  sin  ^  sin  6)  (9) 

Now  since  r  =  |r|  ,  the  phase  distance  from  the  origin  of  the 
point  Pj^  along  the  direction  ^g,  is  by  simplification  of  (9) 

P  (N  r)  =  -  P  r{cos  6  cos  ^  cos((t>  -  tIj^)  -  sin  6  sin  Cl 
c  §  e  g 

(10) 


Hence  the  ground  wave  component  of  the  electric  field  parallel 
to  the  kth  antenna  at  Pj^  is 


E 


rgP 


cos  (4)  -  x]^)  e 


jca  t 


(11) 


4.  Zenneck,  J.  ,  "The  Propagation  of  Plane  Electromagnetic  Waves  over 

a  Flat  Earth  and  Its  Application  to  Wireless  Telegraphy,  "Annalen  der  Physik, 
Vol.  23,  1907,  p,  846 

5.  Feldman,  C.  B.  ,  "The  Optical  Behavior  of  the  Ground  for  Short  Radio 
Waves,  "  Proceedings  of  the  I.  R  E.  ,  Vol,  21,  No,  6,  June  1933,  pp,  764-801 

6.  Piggott,  W.  R.,  "A  Method  of  Determining  the  Polar  Diagrams  of  Long 
Wire  and  Horizontal  Rhombic  Aerials,  "  Dept,  of  Scientific  and  Industrial 
Research,  Radio  Re s earch  Special  Report  No.  16,  London,  H  M  Stationery 
Office,  1948. 

7.  Wait,  James  R.  ,  "Propagation  of  Radio  Waves  over  a  Stratified  Ground,  " 
Geophysics,  Vol.  18,  1953,  pp.  416-422. 

8.  Jordan,  Edward  C.,  Electromagnetic  Waves  and  Radiating  Systems, 

New  York,  Prentice-Hall,  Inc.,  1950. 


where  w  =  p^c  and  ♦j.g  is  an  arbitrary  phase  constant  used  only  to  relate 
Ej.g  to  the  other  electric  field  components.  Substituting  equation  (10)  into 
(ll)ithe  complete  expression  for  the  ground  wave  component  parallel  to  the 
kth  antenna  at  Pj^  is 


)U 


^■gp  ■  ^rg  ■  V  ® 


|t  ^  “  (cos  6  cos  C  cos  (<})g  *  ni^)  -  sin  6  sin 

(12) 


When  only  the  ground  wave  is  present,  equation  (12)  is  sufficient  to  complet 
the  analysis;  however,  when  sky  waves  are  present, the  component  Egp  is 
also  required.  Referring  to  Figure  3  the  sky  wave  ray  which  is  reflected 
at  R  is  delayed  in  phase  by  an  amount  4  at  Pk  where 

e  =  RPi^-QPl,  (13) 


This  simplifies  to 


Therefore. 


I  -  (^o^  _  Pph  cos  Z^l^s  _  4TTh  sin  . 

sin  ij^s  sin  ipg  \q 

the  total  sky  wave  at  Pk  may  be  written 


2poh  sin  i}js 


(14) 


Es  = 


'vs 


1  -  Rv^ 


-jZp-h  sin 


-i<I> 


J'*’vs  , 


+  Ehs  I  1  +  Rhe 


-jZPoh  sin  ij>s 


(15) 


The  phase  of  Eg  at  Pk  relative  to  the  origin  is  Po(Ns  •  T)  similar  to  the 
manner  in  which  equations  (9).  (10)  and  (11)  derived  the  ground  wave  phase 
Po(Ng  •  r).  The  direction  cosines  of  Ng  are 


fg  =  -  cos  4^5  cos  (|)g' 


rrig  -  -  cos  i|jg  sin  (j)g 


both  direct 
and  reflected 


ng  =  +  sin  4^g  (direct) 


(16) 


ng  =  -  sin  ijjg  (reflected) 
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Now  equation  (15)  has  been  written  entirely  in  terms  of  the  incident  field, 
hence  the  third  direction  cosine  with  the  minus  sign  is  not  used.  (If 
equation  (15)  were  written  in  terms  of  the  incident  and  reflected  components, 
the  difference  in  the  third  direction  cosine  would  contribute  the  term  2Por 
sin  4  sin  *j>g,  however,  this  term  was  derived  by  another  method  as  shown 
by  equation  (14).  J  The  coordinates  of  Pj^  are  given  by  equation  (7),  hence 

Po(Ns  •  r)  =  *Por  I  cos  4  cos  cos  4>s  cos 

+  cos  4  cos  sin  <|)g  sin  -  sin  4  sin  i|^g  1  (17) 
Therefore,  the  total  sky  wave  parallel  to  the  antenna  at  Pk  is 


Egp  -  "^IEvsI  sin  cos  ((}>s  *  nk)  i  -  RvC  ^sj  g 

J^hs|  (t  -  — 


+  Ehs  sin«)s  *  nk) 


1 


18) 


where 

Ng  •  r  =  -  r  [cos  4  cos  ijjg  cos  (ct)g  •  rij^)  -  sin  4  sin  i}^gj  (19) 

When  the  ground  wave  and  the  sky  wave  arrive  at  the  same  azimuth, then 

(j)  =  (f)s  =  <j)g.  The  total  electric  field  parallel  to  the  antenna  at  Pk  is  therefore 

the  sum  of  equations  (11)  and  (18) 


Ep  =  ||E,g|  cos  (<(>  -  e'j<Po^g  '  ""  *>-gl 

+  |e;  J  sin+3  cos  (.^  -  m,)  [l  -  e’J'f  ^  *  *vs) 

+  |EkJ  sin(4.-nk)  [l  tR^e'-i^Pohsin+.j^-KPoNs  '  T  *  «hs)|ej“‘ 


(20) 


15 


where 


Ng  •  r  =  -r  (cos  6  cos  C  cos  (<>  -  %)  -  sin  6  sin  Cl 
Ng  •  r  =  -r  [cos  ti^g  cos  C  cos  (4>  -  r|k)  *  *1^8  si*'  C) 


(  Cg  -  1  sin  iji. 

\  ®  <*^^0/ 


w€/ 


COS^  ipt 


(cg-i^)  sin+, 


*  J 


— ^  )  -  cos^  4'< 
W€o  / 


sin 


'h  = 


sin  ij^s  + 


yi'i- 


(21) 


4.  Derivation  of  the  Array  Pattern 

It  may  be  assumed  that  in  an  elemental  length  drk  of  the  kth 
antenna,  at  the  point  Pki  the  parallel  electric  field  Ep  will  induce  a  voltage 
such  that 


dVp  =  Epdrj^ 


(22) 


Thus  the  potential  gradient  is 


drk 


(23) 


The  current  and  voltage  distribution  on  the  Beverage  antenna 
may  be  expressed  by  transmission  line  equations9.  If  the  transmission 
line  (or  Beverage  antenna)  is  immersed  in  a  field  Ep,  parallel  to  the  line, 
the  usual  transmission  line  equation  for  the  potential  gradient  must  include 
a  term  such  as  given  by  equation  (23).  Therefore,  the  transmission  line 
equations  when  written  for  a  Beverage  antenna  are’ 


9.  Beverage,  Harold  H.  ,  Rice,  Chester  W.,  and  Kellogg,  Edward  W., 

"The  Wave  Antenna,  a  New  Type  of  Highly  Directive  Antenna,  "  Transactions 
of  the  American  Institute  of  Electrical  Engineers.  Volume  42,  February  1923, 
pp .  215-266. 
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dV 

drk 


dVn 

(R  +  juL)  I  +  — P 
drk 


^  -  {C  +  juC)  V 

drk 


(24) 


These  equations  may  be  solved  in  the  usual  manner  for  the 
receiving  end  current  as  was  done  by  Herlitz^^,  or  the  line  may  be  assumed 
to  be  terminated  in  the  characteristic  impedance  Zq  ”  ^o  ^  j^O' 
end  (the  point  where  rk  "  rL  for  the  kth  antenna),  and  also  at  the  sending 
end  (rk  =  ro),  and  the  current  at  the  sending  end  found  by  integration  as  was 
done  by  Kellogg^  Either  method  will  yield  the  same  result.  Following 
Kellogg's  procedure,  the  elemental  current  at  Pk  due  to  Ep  is 


dl 


k(-) 


^inl  +  Zin2 


(25) 


where  the  denominator  (Zini  +  Zip2)  is  the  complete  series  impedance  of  the 
circuit.  (Zinl  and  Zin2  are  the  impedances  seen  looking  into  the  line  in  each 
direction;  with  respect  to  the  current  in  the  antenna  conductor,  these 
impedances  are  in  series  through  the  ground.)  Z-j^^  and  are  in  general 

given  by  the  usual  equations 


Zl  cosh  7  (r^  •  r)  +  jZ^  sinh  7  *  r) 

Zq  cosh  7  (rL  *  r)  +  jZl  sinh  7  (rL  r) 


^inZ  "  ^o 


Zj.  cosh  7  (r  -  r^)  +  jZ^  sinh  7  (r  •  r^) 
Zq  cosh  7  (r  •  r^)  +  jZ^.  sinh  7  (r  -  r^) 


25b) 


10.  Herlitz,  Ivar,  "Appendix  B,  Analysis  of  Action  of  Wave  Antenna,  " 

(An  Appendix  to  the  paper  by  Beverage,  Rice  and  Kellogg,  "The  Wave 
Antenna,  a  New  Type  of  Highly  Directive  Antenna"),  Transactions  of  the 
American  Institute  of  Electrical  Engineers,  Volume  42,  February  1923, 

pp.  260-266. 

11.  Beverage,  Harold  H.  ,  Rice,  Chester  W.  ,  and  Kellogg,  Edward  W.  , 

"The  Wave  Antenna,  a  New  Type  of  Highly  Directive  Antenna,  '  Transactions 
of  the  American  Institute  of  Electrical  Engineers,  Volume  42,  February  1923, 
p.  224,  footnote  10. 
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The  elemental  current  at  the  receiving  end  is  retarded  in  phase  and 

changed  in  amplitude  by  the  propagation  constant  of  the  line, hence 


(Note  again  that  rj^  is  the  distance  measured  from  the  point 
(z  =  h,  y  =  0,  X  =  0)  to  Pk  whereas  r  is  the  distance  from  the  origin  to  Pk» 
or  rk  =  r  cos  J  The  total  current  at  the  receiving  end  is 


I  =  -  /  dl  =  r  dl  (27) 

k(-)o  J  k(-)o  J  ^ 

ro 

Substituting  (25)  into  (27)  yields  the  current  at  the  receiving  end 


^inl  ^‘inZ 


(28a) 


where  Ep  may  in  general  be  given  by  equation  (20), 


The  integration  of  equation  (28a)  is  greatly  simplified  if  it  is 


assumed  that 


^■inl  =  ^inZ  = 


(28b) 


This  assumption  is  in  accord  with  the  usual  manner  of  terminating  the 
antenna,  although  in  practice  it  is  never  possible  to  achieve  a  termination 
exactly  equal  to  Zq  for  two  reasons.  First, the  ground  resistance  of  the 
ground  rod  cannot  be  made  zero,  and  second, Zq  has  a  small  reactive 
component  which  cannot  be  matched  over  a  wide  range  of  frequencies. 
Nevertheless,  it  is  a  good  approximation  in  practice,  and,  when  equation  (Z8b) 
is  closely  approximated, equation  (28a)  becomes 


‘k(-)o=H-  f  ■  >^0)  dru 


(28c) 


In  most  instances  it  will  be  true  that  h  «  rQ  so  that  cos  C  ~  1 
and  sin  C  ~  0-  Hence  equation  {28c)  may  be  rewritten  replacing  rk  by  r 


%(-)o 


Epe-7{r  -  ro)dr 
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y  and  the  phase  terms  for  equation  (20)  may  be  rewritten  as 
(Ng  •  r)'  =  -  r  cos  6  cos  (<p  - 

(30) 

(Ng  •  r)'  =  -r  cos  4<s  cos  (({>  -  %) 

Substituting  equations  (30)  into  (20)  yields 

Ip  =||E,g|cos  («  •  <'!>  '  '  *rgl 

tlE„|sin  cos  (6  -  .11,)  [l  -  R„e' '  '^vsl 

+|EKjsin(4.  -  TlJ  f  I  t  '•'s'l  ej'  ^o-- 'I's  (f  '  Ikl  '  ‘hs  I  |eJ“' 

'■  •’  (il)-* 


Hence  the  receiving  end  current  for  this  case  is  given  by 
equation  (29)  with  (31)  substituted  for  Ep.  This  yields 

Ik(.)o  =  if;]  l^rglcos  (4.  -  nk)e-i*rg  ,.ilPor  cos  6  cos  (4  -  nk)  I  -7(r-  ro).j, 

+  [|Evs|sin)|>gCOs(4)- ni^)  (1  -  R^e' sin  e*  j'^vs +|Ehs|sin  (c))  -  qk)(  H 

1  *  r  T 

+  ^s)  e‘j*^hs  •  J  gj^  cos  4^s  cos  ((j)  -  1 -7(  r  *  j. 

To 

(32) 

The  integrals  are  all  of  the  form 


/’ 


L 


.air  +  32  jjj. 


a  I  =  jPoCos  (6,  ijjg)  cos  ((|)  -  nk)  -  7 
a2  =  7ro 


where 
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Therefore 


,  ^jPorLco8(6, cos(<j)  -  nk)  -  ^(rL  *  ro)  jPoi*o  cos (6,  cos (<?  -  nk) 
J  -  e__ _ -e _ 

Tq  jPo^os  (6. 4^5)  cos(4) -Tik)  -  7 

Substituting  the  above  solution  into  equation  (32)  and  factoring  the  common 
terms  of  the  sky  wave  component  yields  for  the  receiving  end  current  of  the 
k(-)th  antenna 


‘k{-)o 


eJwt  r  ,  /  1  .  ejPoLcos6cos((j)  -  nk)  -  7L\  «  1 

=  77-^  E„  cos(4>-nk)e  ^  "S  - 7 - 7 - 77 - — 

2Zo  M  *^8'  ^  \  7  -  iPqCOs  6  cos(<}>  -  Tik)  / 


l^vsl®*” ■ ’^k)  (^  ■  ^v®  ^|Ehs|sin((j)  -  ni^)  ■*■ 


jporocos  »|^s  cos  (<p  -  qj^) 


gjpoLcos^'s  cos(({)  •  nk)  *7L 

7*  jpQCOsi|igCOs((|)  -  Tij^) 


The  current  from  the  kth  pair  will  be 
Iko  “  lk(-)o  +  lk(+)o 

where  lk(-)o  given  by  equation  (33)  and  ik(+)o  given  by  equation  (33)  with  cos 
(<t)  +  0}^)  substituted  for  each  cos(4)  -  Hj^)  term,  and  sin((t)  +rij^)  for  each  sin(<))  -  rj^)  term. 

Introducing  arbitrary  amplitude  and  phase  parameters 
and  Pk(+)  each  antenna  to  correspond  to  an  arbitrary  combining  network 
equation  (34)  becomes 


'k  =  Ck,.,  ly.jo  t  Gk,.,  Ik(^)o  (35) 

In  general  Cfk(.)  /  ^k(+)  Pk(+)  ~  f^k(-)‘  instance, to 

investigate  a  single  antenna  at  =  O.let  =  1,  G^^j  -  0  and  =  0. 

To  investigate  a  difference  pattern  of  two  equal  antennas  let  Gk(-)  =  Gk(+)  and 

Pk(^)  =  0-  Pk(+)  = 
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The  sum  of  a  circular  array  is  therefore 
n 

^total  =  X 
array  k  =  1 

Equation  (36)  is  the  complete  three-dimensional  pattern  of  the 

array. 
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B.  Specification  of  the  Beverage  Antenna  Elements 
1 .  Antenna  Constants  and  the  Wave  Tilt  Angle 

In  order  to  use  equation  (36),  it  is  necessary  to  specify  the 
electric  fields  incident  on  the  array,  the  array  geometry,  the  propagation 
constant  7  of  the  antenna,  and  the  earth  constants.  It  may  also  be  prefer¬ 
able  in  some  cases  to  specify  the  wave  tilt  angle  determined  by  the  earth 
constants.  In  the  case  of  complex  summing  with  phasing,  it  is  also  neces¬ 
sary  to  specify  the  goniometer  amplitude  and  phase  functions  Gj^  and  pj^. 

In  order  to  simplify  these  steps  in  preparation  for  numerical  calculations, 
formulas  for  the  propagation  constant  and  the  wave  tilt  angle  are  required. 

The  propagation  constant  of  a  transmission  line  is  defined  to  be 

7  =  Q  +  Jf3  (37) 

where  the  symbols  are  defined  in  Table  1  and  arc: 

a  =  attenuation  constant  in  nepers/unit  length 

P  =  phase  constant  =  Zir/K  where  \  is  the  wavelength  of  the 
antenna 

The  attenuation  constant  of  the  Beverage  antenna  may  be 
determined  experimentally  and  found  to  vary  with  the  earth  constants  and 
height  above  ground  as  well  as  other  factors.  For  most  purposes  where 
the  Beverage  antenna  is  one  or  two  meters  above  ground,  the  attenuation 
constant  in  the  range  of  3  to  30  mes  will  be  found  to  be  in  the  range  of 
.01  to  .  1  db  per  meter  respectively  for  a  reasonably  good  conductor  such 
as  number  12  or  14  copper  conductor  To  convert  db/meter  to  nepers/ 
meter,  divide  the  db/meter  figures  by  8.686  This  yields  for  the  attenu¬ 
ation  constant  range  mentioned  the  following  very  approximate  but  typical 
values. 


Q  =  10’^  neper  /meter  at  3  me 
a  =  10'^  neper/meter  at  30  me 

These  figures  are  a  good  rule  of  thumb  for  heights  of  about 
one  meter  At  lower  heights  the  attenuation  may  increase  considerably. 

The  phase  constant  may  be  conveniently  expressed  in  terms 
of  a  "velocity  ratio"  n  where: 
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(38) 


An  approximate  value  for  n  may  be  experimentally  determined 
by  exciting  the  antenna  as  an  open  circuited  line,  noting  the  wavelength  on 
the  antenna  and  dividing  by  the  free  space  wavelength.  Then  it  is  possible 
to  write  the  propagation  constant  in  terms  of  measured  data  as: 


-  Q(db/unit  length)  .  i3o( radians /unit  length) 


8,  686 


(39) 


Nearly  exact  values  of  n  or  (3  may  also  be  calculated  from 
Carson^^  and  Wise^3  as  was  stated  previously  for  the  attenuation  constant^'^. 

The  wave  tilt  angle  6  may  be  calculated  from  formulas  given 
by  various  authors^"^  or  it  may  be  measured  experimentally  Since  the 
measurements  are  usually  dependent  on  special  equipment,  calculation  of 
the  tilt  from  the  earth  constants  and  the  frequency  is  the  easier  procedure 
when  a  nonstratified  earth  is  approximated.  Piggot^  gives  the  formula 

1/4 

(esu  units)  (40) 


which  may  be  rewritten  in  MKS  units  as  follows: 


12  Carson,  John  R.  ,  "Wave  Propagation  in  Overhead  Wires  with  Ground 
Return,"  The  Bell  System  Technical  Journal,  Volume  V,  1926,  pp.  539- 
554  (paper  uses  emu  units). 

13.  Wise,  W.  H,  ,  "Propagation  of  High-Frequency  Currents  in  Ground 
Return  Circuits,  "  Proceedings  of  the  IRE,  Vol.  22,  No.  4,  April  1934, 
pp.  522-527  (paper  uses  emu  units). 

14.  Values  of  a  and  n  which  were  obtained  in  a  number  of  experimental 
cases  are  given  in  the  Fifth  Interim  Report  of  this  contract  dated  7  February 
1962,  "Use  of  the  Beverage  Antenna  in  Wide  Aperture  High  Frequency 
Direction  Finding.  " 

15.  Gill,  E  W.  B.  ,  "A  Simple  Method  of  Measuring  Earth  Constants,  " 
Proc.  of  the  Inst,  of  Electrical  Engineers  (London),  Vol.  96,  March 
1949,  pp.  141-145. 
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In  MKS  units  Cq  =  8.  854  X  10^^.  As  a  matter  of  interest 
equation  (41)  is  plotted  in  Figure  4  for  a  variety  of  typical  conditions. 
Typical  values  of<Tg  may  be  inferred  from  these  curves,  other  typical 
values  for  the  earth  constants  are  given  in  Table  2.  It  is  evident  that 
wave  tilt  angles  are  in  general  less  than  20®  with  perhaps  10®  a  rough 
rule  of  thumb.  Note  that  equation  (41)  gives  the  magnitude  of  the  angle 
and  does  not  include  information  relating  the  phase  of  the  tilted  component 
to  the  vertical  component.  Such  data  are  available  in  the  literature®  ^ 

Equations  (39)  and  (41)  may  now  be  used  to  determine  numeri 
cal  values  for  substitution  into  the  pattern  equations. 

2.  Theoretical  Calculation  of  Characteristic  Impedance,  Input 

Impedance  and  Propagation  Constant 

Nearly  exact  values  of  the  characteristic  impedance  and 
the  propagation  constant  7  may  be  calculated  from  the  antenna  geometry 
and  the  earth  constants  by  a  method  devised  by  Carson^ ^  and  modified 
by  Wise.  ^ ^ 


The  exact  values  of  Zq  and  7  for  a  uniform  transmission  line 
are  usually  defined  as^^ 


/  R  +  jeoL 

0  V  G  +  jwC 

(42) 

=  "s/  ( R  +  ju)L)  ( G  +  jojC) 

(43) 

Carson^^  provides  the  following  expression  for  (R  +  jcoL)  in 
electromagnetic  cgs  units  ( abohms/centimeter) 


R  +  jwL  =  rg  +  jZcof  n  f  —  ^  -h  4u)  J  (n/p^  +  j  -  p)e‘2h'pdp 

^  ^  ^  0  (44) 

16.  Skilling,  Hugh  Hildreth,  Transient  Electric  Currents,  McGraw-Hill 
Book  Company,  Inc.,  New  York,  1952. 


TABLE  2  EARTH  CONDUCTIVITY  CONVERSION  TABLE 
(With  Typical  Dielectric  Conatanta) 


Conductivity  (e) 

Reaiativity 

(;) 

Dielectric 

Conatant 

(Typical) 

Type 

of 

Earth 

emu 

(abmho 

cm) 

eau 

(atatmho 

cm) 

MKS 

(mho- 

meter) 

MKS 

(ohm- 

meter) 

(Relative 

Unita) 

Sea  Water 

5  X  10**1 

4  5  X  lolO 

5 

.  2 

81 

Sea  Water 

3  X  10-H 

2.7  X  lolO 

3 

.  33 

81 

Wet  Rich  Soil 

3  X  10-13 

2, 7  X  108 

.  03 

33 

15-36 

Average  Soil 
(Wet) 

1  X  10-13 

9  X  10‘7 

.01 

100 

10  -  25 

Average  Soil 
(Dry) 

3  X  10-14 

2.  7  X  107 

3  X  10-3 

333 

10  -  15 

Poor  Soil 

1  X  10-14 

9  X  108 

1  X  10-3 

10^ 

10 

Poor  Soil 
(Dry) 

3  X  10-15 

2.  7  X  lo6 

3  X  10-4 

3.  3  X  10 

8 

Dry  Sand 

1  X  10-15 

9  X  105 

1  X  10-4 

104 

5 

Dry  Granite 
(Subsurface) 

1  X  10-18 

900 

10-7 

lo'^ 

Probably 

<5 

FIGURE  4 

WAVE  TILT  ANGLE  VS.  FREQUENCY, 
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where  (j.  is  an  integration  parameter  and  has  nothing  to  do  with  perme 
ability. 


The  term  r^  represents  the  distributed  loss  resistance  of  the 
antenna  conductor, and  the  second  term  the  distributed  inductive  reactance  of 
the  antenna  conductor  alone  over  a  perfectly  conducting  ground  plane.  The 
integral  term  formulates  the  effect  of  the  finite  conductivity  of  the  earth. 
Carson  provides  the  solution  to  the  integral  neglecting  displacement  current 
(neglecting  the  dielectric  constant  of  the  earth)  as 

R  +  jwL  =  Tg  +  j2a)f  n  +  4w(P  +  jQ)  emu  (45) 

This  equation  may  be  rewritten  in  MKS  units  throughout  as 

follows 

R  +  jcjL  =  r.  +  j2u,10*‘7f  n  f  +  4w(P'  +  jQ  )  (46) 

®  \  b  /  meter 


whe  re 


P'  t  jQ'  =  (P  +  jQ)10-'^  (47) 

P  and  Q  as  derived  by  Carson  are  complicated  functions  (to 
be  considered  below)  of  rj,  where 

r.  =  2h\/4TT(r  u,  (^^^0  X  10*^)  (48) 

^  s 

In  equation  (48),  when  h  and  (Tg  are  in  MKS  units,  r^,  is  in  emu 
so  that  reference  to  Carson  s  original  equations  for  P  and  Q  may  be  made 
without  conversion.  Equation  (48)  may  be  simplified  by  writing 

r^  ---  .00317  h  (49) 

This  function  is  plotted  in  Figure  5  for  various  conductivities 
in  order  to  illustrate  that  practical  HF  values  for  r^  will  be  in  the  vicinity 
of  .  5  <  Tj,  <  6. 0. 

Wise^^  states  that  above  60  kc  the  dielectric  constant  cannot 
be  neglected.  For  earth  constants  corresponding  to  a  "highly  conducting 
earth,  "  the  error  does  not  appear  to  be  serious  below  10  me.  Above  this 
frequency  the  correction  may  become  significant.  The  effect  of  Wise  s  cor¬ 
rection  is  to  multiply  Carson's  parameter  r^  by  a  complex  correction 
factor  s  where 
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FREQUENCY  IN  MEGACYCLES 
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8  =  t  e 
’w 


ji  J, 


1/2 


emu 


(50) 


The  function  s  is  plotted  in  Figure  6  for  typical  conditions  of  er_  =  .03,  Cj.  = 

12.  It  seems  evident  that  the  correction  should  be  made  at  higher  frequencies 
or  lower  conductivities.  Wise's  correction  obtains  a  new  parameter  rj.^y  = 
rj.8,  which  is  to  be  used  in  Carson's  original  equations  in  place  of  r^..  The 
new  parameter  is  therefore 


r  =  r 
cw  c 


1  +  J 


nl/2 


2c\cr 


g 


emu 


(51) 


This  may  be  rewritten  with  the  ter  m  on  the  right-hand  side  in  MKS  units  as  follows 

Ua  -  1)  1''^ 


^cw  ”  ^c 


1  +  j 


_s_ 

2c\aglO-‘7j 


emu 


(52) 


Equation  (52)  provides  the  same  numerical  value  for  r^,^  as  equation  (51), 
if  MKS  input  data  for  \  and  (Tg  are  used.  This  allows  the  direct  use  of 
Carson's  original  curves  or  equations  without  further  concern  over  con¬ 
version  of  units  (the  papers  by  Carson  and  Wise  are  in  emu  units).  The 
corrected  formula  for  P  +  jQ  is  therefore  given  by  the  following. 


p  +  jQ  =  2(i  -S4)4('''(r 


0-H 


+  -  a,  +  j 


£n 


78Il)(rcJ 

2 


1.7811  r^w 


)(1  -  S4) 


where 


'n  8  ^  I 


1  1 

(  ^cw  ^ 

1'.  '  1 

(  ^cw 

1:2: 

^  2  ) 

3:4' 

^  2  / 

(53) 


(54) 


17.  Equation  (53)  is  suitable  for  self-impedance  calculations  since  Carson's 
parameters  0,  and  are  zero  and  hence  have  been  omitted  As 
Carsonl*^  shows,  it  is  possible  to  make  mutual  impedance  calculations  in 
which  case  these  parameters  are  not  zero.  The  complete  formulas  are 
given  in  the  Appendix  reprint  of  Carson's  paper.  In  general  it  will  be  found 
that  mutual  coupling  as  it  affects  input  impedance  may  be  ignored  when  the 
spacing  between  elements  is  greater  than  the  height  above  ground 


MAGNITUDE  OF  THE  CORRECTION  FACTION  S  FOR  TYPICAL  CONDITIONS 
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(55) 


7 


(56) 


^  [  2  4j  1:2.  V  2  /  L  2  3  4  8j3:4:  \  2  / 


<T-.  = 


r3 

cw 


cw 


.1 1 

cw 


3  ^2  c  ^2  .  c2  _  ^2  Q  i2  c2  .  -72  q2  1  i2 


=  i+iti-il— f— )  -[i+i+i+i+i-— 

2  3  6j2'.  3'.  \  2  /  I  2  3  4  5  10j4'51\  2  / 


—  S4  approximately 


2  / 

+  .  .  . 

(57) 

13 

(58) 

'_L./ 

*  ’"cw  \ 

4'  51  \ 

.  2  / 

(59) 

When  the  dielectric  constant  is  neglected  (Cg  =  1),  the  above 
series  and  equation  (53)  reduce  to  Carson's  original  solution.  Plots  of  P 
and  Q  taken  from  Carson's  original  paper  (in  terms  of  not  r^,^)  are 
given  in  Figures  7  and  8. 

Unfortunately,  Wise  did  not  provide  either  curves  of  P  and  Q 
versus  r^^  or  simple  approximations  valid  in  the  range  of  greatest  interest 
particularly  r^,  <  5.  He  does  provide  formulas  for  r^,  <  ,  25  and  rj,  >  5 
Thus,  in  the  important  range  of  .  25  <  re  <  5  it  is  necessary  to  use  equa¬ 
tion  (53),  or  if  the  correction  by  Wise  can  be  considered  negligible,  then 
the  curves  of  Figures  7  and  8  may  be  used  In  the  latter  case  (negligible 
dielectric)  approximate  formulas  may  be  derived  empirically  for  the  range 
of  .  25  <  r^  <  5.  It  is  then  possible  to  write  a  set  of  simple  equations  for 
all  values  of  r^  as  follows.  (Note  the  range  limits  imposed  by  Carson  have 
been  relaxed  somewhat  since  it  appears  that  his  limits  were  quite  conser¬ 
vative.  ) 


0  vs.  fj 

(TAKEN  FROM  CARSON 'S 
PAPER  FIGURE 3,9^0) 


FIGURE  6 
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P  =  2.54  34  0<rc<6.0  (empirical)  (60) 

.707  r-  -  1 

P  - - 2 -  6.  0  <  r^.  <  00  (Carson)  (61) 


Q  =  -.0368  + 


4.24 


0  <  Tj.  <  .  5  (Carson)  (62) 


Q  = 


_ \ _ 

.  560  +  1  42  r^. 


5  <  r(.  <  6.0  (empirical) 


(63) 


707 

Q  =  - -  6.0  <  Tj.  <  00  (Carson) 


(64) 


where  >  1  5,  a  simple  formula  results  which  may  be  used  with  either  r, 
or  r. 


cw 


P  +  jQ  = -^ -  (I  +  j)  15  <  r^,  <  00  (  65) 

^c 

Up  to  this  point  we  have  been  evaluating  the  term  R  +  juL  in 
equation  (43).  The  term  G  +  jwC  in  equation  (43)  may  be  evaluated  by 
noting  that  in  all  practical  cases  G  0  (nonconducting  support  posts)  so 
that 


G  +  jwC  =  jui 


farads 

meter 


(66) 


Equation  (66)  is  also  an  approximation  as  was  shown  by  Wise 
in  1948.  To  see  that  it  is  just  an  approximation,  one  has  only  to  imagine 
the  earth  turning  into  air  in  which  case  the  height  h  ceases  to  have  signifi¬ 
cance.  Nevertheless,  the  approximation  is  a  good  one  for  typical  conductiv 
ities,  and  no  analysis  of  a  more  exact  solution  will  be  provided.  Wise’s 
paper  is  reproduced  in  Appendix  IV,  should  further  reference  be  required. 

18.  Wise,  W.  H.  ,  "Potential  Coefficients  for  Ground  Return  Circuits,  " 

The  Bell  System  Technical  Journal,  Vol.  27,  April  1948,  pp.  365-371. 
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The  quantity  rg  in  equation  (46)  may  be  determined  from  any 
suitable  text.  For  instance 

ZKn/T  ohms 

r-  =  — r -  — : — ; - -  (67 

®  b  unit  length 

where  K  =  41.6  X  10"^  for  copper^^ 

Combining  the  solutions  found  thus  far  in  equations  (46),  (66) 
and  (67)  the  characteristic  impedance  is  given  by  equation  (42)  as 


_  .832  ~+j2ufn(^)+  4u.(P  +  jQ) 

juZ  ttCq 


(MKS) 


The  propagation  constant  is  given  by  equation  (43)  as 


/  n/I  /2h\  \  10' 


(MKS) 


The  input  impedance  of  the  Beverage  antenna  may  also  be 
calculated  from  the  above  relations  for  Zq  and  7  using  the  familiar  trans¬ 
mission  line  equation 


Zt  cosh  7L  +  jZ  sinh  7L 

2.  =2  — — - - - 

°  Zq  cosh  7L  +  jZ^  sinh  7L 


Note  that  to  obtain  an  ideal  (reflectionless)  termination,  the 
terminating  impedance  must  have  a  reactive  component.  This  probably 
partially  accounts  for  the  fact  that, in  practice,  it  has  been  impossible  to 
find  a  value  for  the  terminating  resistance  which  eliminates  oscillatory 
variations  in  the  input  impedance  with  frequency  , 

19.  Skilling,  Hugh  H.  ,  Electric  Transmission  Lines,  McGraw-Hill  Book 
Co.  ,  New  York,  1951  Equations  (12-15)  through  (12-17).  K  =  40.2  X 
10"9  for  silver,  83  X  10"9  for  brass,  etc. 


I 
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3.  Estimate  of  the  Radiation  Resistance 


The  question  is  sometimes  raised  as  to  the  radiation  resis 
tance  of  the  Beverage  antenna.  An  approximate  value  may  be  found  by 
reference  to  formulas^®  for  the  radiation  resistance  of  a  matched  two- 
wire  transmission  line  section.  The  following  formula  may  be  used  for 
the  latter 


Rr  =  60(2iJoh)^  . 


cosh  aL 


sin  {2(3qL) 


-ciL 


(71) 


For  a  single  wire  over  a  good  ground  plane,  the  radiation 
resistance  will  be  one-half  the  above  value.  Simplifying  yields 


Ry  =  4730 


cosh  qL 


sin  {2poL) 


.  -  aL 


(72) 


It  is  evident  that  even  the  best  line  (a  =  0)  will  have  a  maxi¬ 
mum  radiation  resistance  when 


4ttL 


wheremis  any  integer. 


This  reduces  equation  (72)  to 

Rj.  =  4730  ohms  (73) 

For  h  =  1  meter  and  \q  ranging  from  100  to  10  meters  (the 
high  frequency  range,  3  to  30  mcs)  the  radiation  resistance  ranges  between 
,47  and  47  ohms.  Thus  the  radiation  resistance  never  exceeds  10  percent 
of  the  characteristic  impedance,  and  is  therefore  probably  not  an  important 
factor  in  determining  the  input  impedance  of  low  Beverage  antennas.  It  is 
also  apparent  that  the  efficiency  of  the  Beverage  antenna  will  be  low  since 
Rj.  is  low  compared  to  Zq. 


20  Storer,  James  E..  and  King,  Ronold  "Radiation  Resistance  of  a  Two- 
Wire  Line,  "  Proceedmg_s_of_tb£j[R^  Vo  1.  39  No.  6,  1951,  pp.  1408- 
1412.  Equation  (5),  Case  I,  first  unnumbered  equation 
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C,  Effective  Height,  Sensitivity  and  Signal-to-Noise  Ratio 


The  effective  height  of  the  Beverage  array  may  be  defined  as^^: 

I  I 

(  ^otal  } 


where  (t  =  0,  |E^g|  =  0,  |  E^g|  =  0,  and  all  Gj^  =  1,  all  Pj^  =  0.  Substituting 
equations ( 33)  and (36)  into  (74)  yields  for  the  above  special  case 


^  I 

I  vg|  k=l 


1  .  e(jPoCOs6  cosTik  -  7)L 
T  •  Pq  cos  6cos  Hk 


It  is  evident  from  Figure  2  that 


Erg  I  =|Evgl  tan  6 


Hence  the  effective  height  is 


V  tan  6 

he-  — 


1  -  e(jPo  cos  6  cos  nk  -  7) 
Y  -  Pq  cos  6  cos 


This  may  be  rearranged  as  follows 


Ltan  6 


Z  cos,, 

k=  1 


(^) 


jPo  cos  6  cos  r|k  -  7)  L/ 2 


where 


^k  =  7  -  jPoCOs6  cos  r|j^  (7< 

Since  the  E/  2)  |  ^  unity,  equation  (  7  8)  simplifies  tc 


21.  Schelkunoff,  Sergei  A  ,  and  Friis  Harald  T.,  Antennas  Theory 
and  Practice,  New  York,  John  Wiley  &  Sons,  Inc  ,  1952,  p.  301,  problem 
9,13.  It  is  true  that  Eyg  is  not  parallel  to  the  antenna  conductors;  however, 
it  is  more  convenient  to  refer  the  effective  height  to  a  vertically  polarized 
field  normal  to  the  earth's  surface 
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k=  1 


cos  • 


-7L/2 

e 


Since  7  =  a  +  jP, equation  (80)  simplifies  to 


hp  =  —  tan  6 
2 


aL/2 


Z 

k=  1 


(80) 


(81) 


If  equations  (37)  and  (38)  are  substituted  into  equation  (79)  and  the  result 
simplified  and  substituted  in  equation  (81),  and  also  a  =  0  the  effective 
height  becomes 


he 

Q  =  0 


L 

2 


tan  6 


Z 

k=  1 


cos  r)k 


cos  6  cos 


cos  6  cos 


(82) 


For  a  single  antenna  atiQj^  =  O.this  becomes 


he 
a  =  0 
single 
antenna 


Y  tan  6 


sin 


Lit  /  1 


^o 


—  cos  6 
n 


a"! 

r —  1  —  cos  0 
_^0 


(83) 


When  the  velocity  ratio  n  =  1  and  cos  6  approaches  unity,  then  (sin  X)/X  "  1 
and 


he  = 


tan  6 


(84) 


This  may  be  taken  as  a  rule  of  thumb  for  the  effective  height  of 
an  ideal  Beverage  antenna. 


It  will  be  noted  however  that  when  cos  6  approaches  unity  6 
approaches  zero,  hence  tan  6  approaches  zero.  Hence  equation  (84)  is 
only  valid  for  small  wave  tilt  angles,  and  since  tan  6  =  6  when  6  is  small 
equation  (84)  may  be  written  as. 
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2 


(85) 


where  is  in  radians. 


Equation  (83)  implies  a  minimum  effective  height  at  some  frequency. 
When  (81)  is  simplified  to  a  single  antenna  then 


u  —  *  ,  -qL/2 

hg  -  Y  tan  0  e 

Q  s  0 

single 

antenna 


sinhfr- jPocos6)  2] 


(Y*  jPoCos  ^)  ^ 


(86) 


According  to  equation  (83)  a  minimum  effective  height  for  a  =  0  occurs 
when 


mTT 


where  m  may  be  any  integer.  Thus  the  frequency  of  the  first  minimum 
hg  is  approximately 


f  = 


nc 


(i.cos6) 


L(  1  -  n  cos  6) 


(87) 


or 


f  = 


300n 


L  ( 1  -  n  cos  6) 


megacycles 


(88) 


This  may  be  thought  of  as  a  cutoff  frequency  for  the  antenna,  for  only 
below  this  frequency  can  relatively  large  effective  heights  be  obtained 
In  order  to  design  to  maintain  this  frequency  above  30  mcs,  it  is  neces 
sary  that  f  >  30  and  hence 


10  n  =  L  ( 1  -  n  cos  6)  =  L  -  nL  cos  6 
n(  10  +  L  cos  6)  =  L 

-  L 

10  +  L  cos  fi 


(89) 
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but  n  >  I  cannot  be  obtained  without  modifying  the  antenna  so  the  length  is 
limited  to 

L  =  ( 10  +  L  cos  f')  n 
M  i  -  n  cos  ^)  =  lOn 


L  = 


lOn 


for  n  <  I 


1  -  n  cos 

To  successfully  use  300*meter  antennas, then 


30  = 


1  -  n  cos  6 


or 


30  -  30n  cos  6  =  n 


cos  6  = 


30  -  n 
30n 


30  29 

which  means  for  ^  =  0,  30n  =  30  -  n,  n  =  —  (min  n)  or  for  n  =  1  cos  ^  =  — 

Therefore 


^max  ■  max  for  n  =  I 


n  =  .  965  min 


(90) 


Experience  has  shown  for  h  =  I  meter  these  criteria  can  be  met  for  L  -  300 
at  f  =  30  mcs. 

The  sensitivity  of  the  Beverage  antenna  may  be  estimated  from  the 
minimum  field  strength  which  may  be  received,  determined  by  noise  limita¬ 
tions.  For  this  estimate  we  will  make  the  simple  assumption  that  the 
antenna  equivalent  circuit  consists  of  the  radiation  resistance  in  series 
with  a  loss  resistance  equal  to  the  characteristic  impedance  of  the  antenna. 
The  noise  arising  in  the  antenna  may  be  calculated  by  Nyquist's  theorem 
which  states  that  the  mean  square  fluctuation  e.m.f.  in  the  frequency  band 
V  to  V  +  dv  appearing  in  an  impedance  Z  =  R  +  jX  is  given  by'^'^ 


22.  Nyquist,  H. ,  "Thermal  Agitation  of  Electric  Charge  in  Conductors,  " 
Physical  Review,  Vol.  32,  July  1928,  p.  110. 
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de  =  */4WRdv  (91) 

where  k  =  Boltzmann's  constant  =  1. 374  X  10"^^  joules  per  degree  Kelvin, 
and  T  is  the  absolute  temperature  in  degrees  Kelvin.  The  value  R  is 
assumed  to  be  constant  over  the  bandwidth  dv;  if  R  is  a  function  of  frequency, 
a  more  general  form  of  equation  (91)  is  required.  However,  for  the 
present  analysis  we  will  be  concerned  with  bandwidths  which  are  less  than 
one-half  of  one  percent  of  the  center  frequency,  so  that  while  the  previously 
calculated  radiation  resistances  are  rapidly  changing  with  frequency  the 
changes  are  small  over  the  bandwidths  considered  and  equation  (91)  is 
closely  approximated 

Burgess^*^  shows  that  the  radiation  resistance  of  an  antenna  can 
apparently  be  the  source  of  thermal  noise  and  that  it  obeys  Nyquist's 
theorem  when  the  aerial  is  in  radiative  equilibrium  with  its  surroundings. 
Thus 


der  =  N/4kTRj.dv  (92) 

The  equivalent  noise  temperature  Tj.  of  the  radiation  resistance 
Rj.  is  defined  as  that  temperature  which  must  be  ascribed  to  Rj.  so  as  to 
give  the  value  of  the  received  noise  electromotive  force  when  substituted 
in  Nyquist's  equation,  so  that 

=  ^y4kTJ.Rrd  (93) 

Thus  an  antenna  in  an  enclosure  at  a  uniform  temperature,  T,  is  in 
effect  the  source  of  a  thermal  e  m.f.  whose  value  corresponds  to  the 
condition  Tj.  =  T.  If  the  antenna  temperature  T-^  is  not  equal  to  T  of  the 
enclosure,  radiative  equilibrium  does  not  exist  and  there  will  be  an 
unbalanced  energy  flow  between  the  antenna  and  its  surroundings  in  the 
sense  which  tends  to  equalize  T^  and  T.  It  is  concluded  that  the  radiation 
resistance  of  an  antenna  in  free  space  has  an  equivalent  noise  temperature 
of  zero;  expressed  otherwise, there  is  no  other  source  of  radiation  present 
to  return  energy  to  the  antenna  and  thus  to  induce  a  fluctuation  e.m  f 

23  For  instance  see  Terman,  F.  ,  "Radio  Engineers  Handbook,"  p,  476. 

24.  Burgess,  R  E.,  "Noise  in  Receiving  Aerial  Systems  "  Proceedings 
of  Physical  Society,  VoL  53,  May  1941,  p.  293. 

25.  Burgess  points  out  that  this  statement  is  in  disagreement  with  other 
authors  in  papers  published  prior  to  his  1941  paper.  The  validity  of 
Burgess'  paper  now  seems  well  established  however. 
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Ty  will  in  general  be  independent  of  the  bandwidth,  but  it  will  not 
be  independent  of  antenna  gain  or  effective  height  if  the  noise  source  is 
not  uniformly  distributed  in  all  directions.  Burgess  points  out  therefore 
that  Tj.  will,  in  general,  exceed  T,  and,  in  fact,  high  frequency  measure¬ 
ments  by  Jansky  have  shown  that  Tj.  is  always  greater  than  T  and  can 
exceed  lOOOT  in  unfavorable  circumstances. 

Assuming  that  the  proper  Tj.  is  attributed  to  the  radiation  resistance 
then  the  noise  e.m.f.  of  the  lossless  antenna  is 

ej.  =  V  4kT  (94) 

where  B  is  the  bandwidth.  The  noise  e.m.f.  of  the  assumed  Beverage 
antenna  is  therefore 

er  =  ^^4kfrB(■RI.  +  Rq)  (95) 

Now  the  voltage  output  of  an  antenna  is  related  to  an  incident  signal  field 
E  by 


V  =  Eh, 


(96) 


Hence  for  unity  signal -to-noise  ratio  at  the  output  terminals,  the  incident 
signal  field  strength  must  be 


E  = 


N/4kTr(Rr  +  Rq) 

he 


(97) 


Thus  for  typical  values  of  Rj.  =  lOS?,  Rq  =  400P,  hg  =  10  meters  B  =  3  X 
103  cycles  and  T  =  300  degrees  Kelvin, one  could  expect  the  unity  signal-to- 
noise  ratio  to  occur  at 

E  =  n/(4)(1.374  X  10-23)(300)(41o)(3  x  103)  pv 

10  meter 


or  E  =  .014  microvolt  per  meter.  Of  course,  the  "signal"  in  this  case 
could  be  atmospheric  noise,  or  any  other  external  source. 
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D.  Beverage  Antenna  Array  Calculations  with  a  Digital  Computer 

To  differentiate  the  present  calculation  from  past  calculations  of  a 
similar  nature  on  this  and  previous  projects,  the  present  computer  pro¬ 
gram  is  referred  to  as  the  "Beverage  Array  11^.  Program."  The  program 
is  written  for  a  GE-<i25  using  IBM  cards  and  the  GE  "WIZ"  program 
language.  (Note:  Much  of  the  information  which  follows  in  this  section  is 
of  interest  only  to  personnel  at  this  laboratory.  It  is  provided  herein  in 
order  to  maintain  a  permanent  record  and  for  the  future  convenience  of 
project  personnel.  Those  readers  not  interested  in  the  computer  program 
ming  details  may  skip  to  Section  E  for  plotted  patterns.) 


The  Beverage  Array  II^.  Program  calculates  the  response  of  a 
Beverage  antenna  or  a  circular  array  of  radial  Beverage  antennas  to  the 
following  incident  waves. 


(1) 

(2) 

(3) 


Vertically  polarized  ground  wave 
Vertical  polarized  sky  wave  V 

Horizontally  polarized  sky  wave 


At  the  same  azimuth 
if  more  than  one  is 
specified  at  one  time 


The  equation  calculated  gives  the  current  at  the  receiver  end  of  a  Beverage 
antenna  or  the  array  under  excitation  by  any  or  all  of  the  above  waves  The 
equation  calculated  is  equivalent  to  equation  (36)  and  is  written 


N  =  10 
K  =  4N  +  1 


=  I 


GKeJP(K  t  1) 


1  ^'^2 


(T,  +  T,)T 


i  • 


(99) 


K  =  0 
N  =  0 


where 


/,  jPqL  cos  6  cos  (4)  -  r](N)]- ^  L 
T2  =  Ej.g  cos  t4>  -  r](N)le  ^8  '  - - 


7  -  jPo  cos  6  cos  [(j)  -  r](N)l 


JPqRq  cos  a  cos  [<t>  -  r](N)) 


— 

sin  i|j  cos  (cf)  -  r](N)]  e  1 


•r  sin 
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^3  = 


sin  Ip 


ll 


,h  sin  Ip, 


JPqL  cos  Ip  cos  (<p  -  n(N)l  ■  *>  L* 
7  -  jPo  cos  ip  cos  (<j)  -r)(N)l 


.  gjPo^o  ‘^°®  4>  cos  {4>  -  n(N)l 

(100) 


Auxiliary  calculations  in  order  to  obtain  the  antenna  current  given 
above  include  the  Carson  impedance  parameters  P  and  Q,  and  the  propaga¬ 
tion  constant  for  the  Beverage  antenna,  see  equations  (53)  and  (69) . 

The  array  of  Beverage  antennas  is  assumed  to  be  circular  with  pre¬ 
scribed  inner  and  outer  radii  (Rq,  R^)  specifying  antennas  of  desired  lengths 
(L).  If  a  phased  array  is  required,  provision  is  made  for  the  multiplication 
of  each  antenna  output  by  a  magnitude  and  phase  number  which  has  the 
effect  of  modifying  the  effective  array  geometry  as  desired  (as  for  instance, 
the  insertion  of  delay  lines).  This  array  of  multipliers  is  referred  to  as 
the  goniometer  function  (  see  equation  (35)]  . 

The  calculation  proceeds  by  summing  antennas  by  pairs,  symmetric 
about  0®  azimuth.  A  center  antenna  at  0  degrees  may  also  be  specified  and 
calculated.  The  antennas  do  not  need  to  be  equally  spaced 

The  array  consists  of  (2N  +  1)  antennas,  or  less,  divided  into  (N) 
pairs  either  with  or  without  the  degenerate  pair  N  =  0,  which  is  a  single 
antenna  located  at  0  degrees.  The  specification  of  the  number  of  pairs  is 
carried  out  by  specifying  the  variable  NIM  (either  0  or  1)  naming  the 
minimum  number  of  the  N  array.  Thus  when  NIM  =  0,  the  single  0® 
antenna  is  included:  when  NIM  =  1,  the  calculation  begins  with  the  first 
true  pair  and  the  center  antenna  is  omitted.  The  total  number  of  pairs 
computed  is,  therefore,  NMAX  -  NIM  including  the  degenerate  center 
antenna.  (This  means  it  is  possible  to  use  NIM  values  greater  than  one, 
however,  there  seems  to  be  no  advantage  to  this.)  The  variable  NMAX 
specifies  the  number  of  antenna  pairs.  The  azimuth  of  each  pair  is  speci¬ 
fied  by  the  variable  ETA(N)  symmetric  about  the  0  degree  axis  There¬ 
fore,  ETA(N)  where  N  =  0  specifies  the  0  degree  antenna 

Although  the  antennas  are  specified  in  pairs,  a  suitable  goniometer 
function  with  0  magnitude  for  the  appropriate  Kth  antenna  will  specify  an 
array  of  irregular  antenna  spacings  as  desired. 


I 

I 
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The  Beverage  array  calculated  by  the  program  is  diagrammed 
schematically  below.  One  antenna,  N  =  Z,  on  one  side  has  been  eliminated 
by  specifying  a  4ero  goniometer  function  at  that  position  on  the  right 
side  of  the  array. 

0* 

K  =  0.  1 
K  =  4,5  N  =  0 


Example  Goniometer  Function 
^(K)< 


N 

0 

1 


lej° 

I.JO 

IcJO 

Oej® 

leJ“ 


K 


FIGURE  9  SKETCH  OF  BEVERAGE  ARRAY 
CALCULATION  GEOMETRY 


Table  3  summarizes  the  computer  input  parameters  and  their 
definitions.  Table  4  gives  the  order  of  input  cards  to  the  program. 
Table  5  gives  the  description  of  sense  switch  options  available  to  the 
program  operator.  Table  6  gives  the  available  output  quantities. 

The  printed  output  from  the  computer  consists  of: 
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TABLE  3 

ANALYSIS  AND  MACHINE  SYMBOLS 
INPUT  PARAMETERS 


Analysis 

Machine 

b 

B 

h 

H 

"g 

Ro 

SIG 

EPG 

RO 

R 

RL 

LAMBDA 

n 

ETA 

Gk(±) 

G(K)  Goniometer 

Pk(±) 

G(K  +  1)  Goniometer 

Erg 

ERG 

rg 

^vs 

FRG 

EVS 

^vs 

FVS 

Ehs 

EHS 

^s 

FHS 

4) 

PHI 

4^ 

PSI 

a4) 

DELPHI 

4Vnax 

PHIMAX 

A 

DELPSI 

4^max 

PSIMAX 

‘I’o 

PHIO 

4^0 

PSIO 

N 

N 

P 

P 

Q 

Q 

a 

ALPHA 

P 

BETA 

B 

BMWTH 

A 

Amperit 

Zo 

ZO 

NOF 

Remark 


Conductor  radius 
Height  above  ground 
Ground  conductivity 
Ground  dielectric  constant 
Array  inner  radius 
Array  outer  radius 
Wavelength 
Location  of  antenna 
Amplitude  of  Antenna 
Phase  of  Antenna 
Amplitude  ground  wave 
Phase  ground  wave 
Amplitude  vertical  sky  v^ave 
Phase  vertical  sky  wave 
Amplitude  horizontal  sky  wave 
Phase  horizontal  sky  wave 
Azimuth  variable 
Polar  variable 
Azimuth  increment 
Maximum  azimuth 
Polar  angle  increment 
Polar  angle  maximum 
Initial  <j) 

Initial  ijj 

No.  antenna  pairs  in  array 
Antenna  parameters  real  part 
Antenna  parameters  imaginary  part 
Attenuation  constant 
Phase  constant 
Beamwidth  criterion 
Amplitude  criterion  for 
Characteristic  impedance 
Number  of  frequencies  to  be  run 
in  parameter  study 


All  input  numbers  in  MKS  units 
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TABLE  4.  LIST  OF  INPUT  CARDS  FOR  BEVERAGE  ARRAY 
RUN  ON  CE-225  WITH  WIZ  LANGUAGE 


Ou«ntity  of 
Dou  on  Numb«r«  on 

Card  No.  Each  Card  Each  Card 

i  T.NOF  Z 

la  o.VR.Rq.Xq  4 

^  B.H.F  J 


4  Evf.^vg.  6 

Eyg* 

^hi'^hi 

5  NIM.NMAX  *  6 


6  MN)  (NMAX^Uwhan 

NIM  ■  0  and  (NMAX) 
when  NIM  ■  1 


Deacrtpiion  of  Card 

Mode  of  caicuiation,  T  •  0  or  T  •  9  faae  opttona  below). 
Number  of  fre^uan^lee  to  be  run  in  parametar  itudy 
when  T  •  5.  uhen  T  •  0  eat  NOF  •  0. 

Not  required  unleii  aanae  ewttch  S  it  down.  When 
raquirad  put  a  in  naperi  per  mater,  vaioctty  ratto 
fdimenaionieii).  R^  in  ohmi  and  in  ohma  where 
Rq  *  j^o  *  ^o  *  charade rlatic  impedance. 

Wire  radiua  in  metera,  wire  height  in  metera  and  fra* 
quancy  in  mca. 

Ground  conatanta.  Conductivity  a^  in  MKS  unita.  and 
fg  lha  dieiectrir  conatant  fdimanaionieaa). 

Wava  parametara.  Vertical  ground  wave,  vertical  aky 
wava  and  horiaontal  aky  wave  in  arbitrary  amplitude 
unita  (E)  and  degraaa  (4). 

Array  parametara.  Whan  cantar  antenna  fthe  dagenar- 
ate  pair)  ia  prcaant  aet  NIM  •  0.  whan  canter  antanna 
la  abaant  aat  NIM  •  1.  NMAX  ia  numbar  of  aymmetricai 
paira,  R^  ia  tha  innar  radiua  In  metara,  R^  ia  tha  outer 
radiua  in  metera.  4^  ia  tha  normaliaad  aaimuth  in 
degraaa  and  ia  tha  normaiiead  eievation  angle  In 
degraea.  (When  calculating  aaimuth  patterna  at  diffarent 
eiavation  angiaa  epecify  the  eievation  angle  with  4n*) 

Locationa  of  paira  of  antennae  In  aaimuth  in  degreea. 
(For  a  pair  of  antennae  at  a4*  put  only  a  eingle  4  on  card, 
for  two  pair  at  a  2*  and  *  4*  put  only  a  2  and  a  4  on  card, 
etc.  When  NIM  ■  0  firet  number  on  card  will  be  aero.) 
Maximum  of  10  paira  piua  the  center  antenna  (or  21 
antennaa). 


6a 


7 


7a 


L 

0 


G(K)  (4NMAX  4  2) 


6 


BMWTH,  2 

AMPCRIT 


Deacription  of  Option 

Random  data  changea  for  aucceaaive 
caeea  (when  T  >  0.  aet  NOF  ^  O). 

Not  uaed. 

Not  uaed,  I  Eliminated  from 
Not  uaed.  j  earlier  program 
Not  uaed.  J 

Only  the  frequency  ia  changed  in 
aucceaaive  caaea  and/or  4n> 

Provide  aeven  carda  aa  uaual  for  the 
firat  caae.  then  the  eighth  card  will 
have  f.  4f|>  '^n 

etc. .  and  thereafter  only  one  card 
per  caae. 


Not  required  unlaaa  aenae  awitch  1  ia  down.  When  aenae 
awitch  i  ia  down  apecify  an  ampiituda  and  phaae  conatant 
for  the  center  antenna  and  ail  other  paira  uaed.  There 
ahouid  be  two  numbera  (which  are  the  firat  numbera  on 
the  card)  to  apecify  the  amplitude  and  phaae  (in  that  order) 
for  the  center  antenna.  Thereafter  there  ahouid  be  four 
numbera  for  each  antenna  pair,  (n  the  firat  group  of  four 
numbera  the  firat  two  numbera  correapond  to  the  right* 
hand  eiement  (aa  aeen  from  the  array  center)  of  the  pair 
cioaeat  to  the  center.  The  aecond  two  numbera  corre¬ 
apond  to  the  left  hand-element  of  the  aame  pair.  The 
firat  number  in  each  group  ia  an  arbitrary  amplitude 
conatant  and  the  aecond  number  ia  a  phaae  conatant  in 
degreea.  The  next  cioaeat  pair  of  antennaa  ia  apecified 
by  another  group  of  four  numbera  and  ao  on.  Note  that 
two  numbera  muat  be  uaed  to  apecify  the  center  antenna 
even  when  it  ia  omitted  by  aetting  NiM  ■  1  on  card  6.  in 
thia  caae.  the  firat  two  numbera  on  card  6a  will  alwaya 
be  00. 

Calculation  parametera.  Initiai  azimuth  point,  azimuth 
increment,  maximum  azimuth  point,  initial  elevation 
point,  elevation  increment,  maximum  elevation  point  (up 
to  180*)<  All  unita  in  degreea. 

Not  required  unleaa  aenae  awitch  3  ia  down.  When 
required  put  beamwidth  (degreea)  and  amplitude  cri¬ 
terion  (normalized  amplitude)  on  thia  card.  For 
inatance,  for  20*  beamwidth  at  the  3  db  pointa  the  2 
numbera  would  be:  (20,  .707). 


TABLE  5 


SENSE  SWITCH  OPTIONS- BEVERAGE  ARRAY  II^. 


Sense 

Switch 


None 


1. 


Z. 

3. 

4. 

5. 

6. 


7. 

8. 

9. 

10. 


_ Resulting  Output _ 

Pattern  calculation  giving  the  following  output  arranged  as 
shown: 


SIC 

EPC 

EVC 

FRG 

EVS 

FVS 

EH3 

FHS 

Antenna  Parameters 

NIM 

NMAX 

^^0 

B 

H 

ETA(N)  ETA(N  +  1) 

‘l>o 

‘l^o 

(j> 

^max 

^max 

Normalized  Point 


IMAG  HE 

PHI  IMAG  IPHASE  INORM  -Ndb  HE 

Goniometer  function  incorporated.  (See  Fig  11)  Note  the  G(K) 
array  is  printed  in  an  array  above  all  the  other  output  data  only 
when  and  always  when  switch  1  is  down.  (Requires  card  6a.) 
Print  a,  p,  VR,  Rq,  Xq  in  addition  to  usual  data 
Beamwidth  option.  (Requires  card  7a  ) 

Print  R^,  R^  Phase,  Rp|,  Rpj  Phase 
Rd  a,  VR,  Rq,  Xq,  input  data  (Requires  card  la  ) 
Self-normalization  (Requires  no  additional  cards.)  Computer 
searchs  in  azimuth  or  elevation  until  first  maximum  is  found, 
pattern  is  normalized  to  this  maximum.  When  function  is 
decreasing  at  zero  the  pattern  is  normalized  at  zero. 

Exit  to  next  case  after  calculating  normalized  point. 

Not  used. 

Unconditional  exit  to  next  case, 

Diagnostic  print. 


TABLE  6 


LIST  OF  PROGRAM  OUTPUT  QUANTITIES 


Analysis 

Machine 

Remark 

a 

ALPHA 

Attenuation  constant 

R 

• 

BETA 

Phase  constant 

n 

VR 

Velocity  ratio  v/c 

R 

R 

Antenna  resistance 

X 

X 

Antenna  reactance 

6  . 

DELTA 

Wave  tilt  angle 

hxl 

IMAG 

Magnitude  of  antenna  current 

‘^T 

IPHASE 

Phase  of  antenna  current 

INORM 

Normalized  magnitude 

-Ndb 

Normalized  magnitude  (db) 

HE(M) 

Effective  height  (meters) 

Rv 

RV(0) 

Magnitude  vertical  Fresnel 
reflection  coefficient 

Pv 

RV(1) 

Phase  vertical  Fresnel 

reflection  coefficient 

Rh 

RH(0) 

Magnitude  horizontal 
reflection  coefficient 

Ph 

RH(1) 

Phase  horizontal  reflection 
coefficient 
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(1)  The  specified  input  parameters. 

(2)  The  calculated  characteristic  impedance,  propagation  constant 
and  wave  tilt  angle,  if  desired, by  placing  sense  switch  2  down. 

(3)  The  Fresnel  reflection  coefficients,  if  desired, by  placing 
sense  switch  4  down. 

(4)  The  initial  current  value  at  the  normalized  point  on  the 
pattern  (both  azimuth  and  elevation)  and  the  effective  height. 

(5)  The  antenna  pattern. 

The  output  numbers  associated  with  each  azimuth  or  angle  of  inclination 
obtained  for  the  antenna  pattern  include  the  magnitude  of  the  current,  the 
phase  of  the  current,  the  current  normalized  to  the  initial  normalizing 
point,  the  normalized  current  in  terms  of  db  down  with  respect  to  the 
normalized  point,  and  the  effective  height  of  the  antenna  or  antenna  array 
(calculated  only  for  the  ground  wave  case).  Both  azimuth  plane  and  polar 
plane  plots  are  available.  An  input-output  block  diagram  of  the  calculation 
is  included  in  Figure  10.  A  detailed  block  diagram  of  the  program  is 
given  in  Figures  11  thru  18. 

One  option  (sense  switch  3)  enables  a  series  of  cases  to  be 
calculated  in  terms  of  beamwidth  and  an  amplitude  criterion.  With  switch 
(3)  down  and  the  desired  half  beamwidth  for  a  given  normalized  amplitude 
specified,  successive  array  patterns  will  be  calculated  but  not  printed  until 
the  amplitude  and  beamwidth  criteria  are  met.  Then  with  switch  (3)  reset 
to  zero, the  desired  pattern  will  be  calculated  and  printed.  For  the  rejected 
cases  the  description  of  the  array  and  the  beamwidth  at  the  prescribed 
amplitude  criterion  is  printed  before  exit  to  the  next  case. 

With  switch  (6)  down, the  program  performs  a  self-normalizing 
calculation  as  follows: 

(!)  Calculate  the  antenna  response  at  4)jyf, 

(2)  Increment  the  appropriate  variable  (4)  for  azimuth,  i|j  for  polar 
plane  or  elevation). 

(3)  Compare  this  contribution  with  the  previously  calculated  point 
until  1(4)  +  ,:^4))  -  I(4>)  is  negative. 


FIGURE  10. 

BEVERAGE  ARRAY  EC  SIMPLIFIED  INPUT-OUTPUT-SENSE  SWITCH  BLOCK  DIAGRAM. 
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The  normalized  point  is  !(<)>).  This  is  a  satisfactory  method  of 
normalization  when  there  is  a  single  main  lobe  adjacent  or  coincident 
with  the  initial 

With  switch  (7)  down.the  normalized  point  is  obtained  only,  in  addi¬ 
tion  to  the  initial  parameters  as  desired  which  are  available  from  switch 
(Z)  and  switch  (4).  This  provides  the  effective  height  and  location  of  the 
main  lobe,  but  exit  to  the  next  case  is  obtained  automatically  so  that  the 
complete  pattern  is  not  calculated.  This  option  is  useful  for  instance  in 
computing  polarization  error  vs  frequency. 

Further  notes  on  program  operation  include  the  following. 

(1)  For  goniometer  function  operation  the  complete  array  beginning 
with  the  0*  antenna  is  called  for,  even  though  G(o)  may  be 
zero.  Therefore,  NIM  =  0  for  all  goniometer  cases  and  r}{o) 

=  0  must  begin  the  ETA(N)  array. 

{Z)  For  operation  without  the  goniometer  function  and  without  the 
center  antenna, it  is  necessary  to  begin  with  the  first  pair  in 
specifying  antenna  locations. 

(3)  For  parameter  studies  with  the  frequency  as  the  only  variable, 

the  input  for  successive  cases  consists  only  of  the  desired  (me) 
frequencies  and  the  desired  normalized  point  (j)pq,  For  these 

cases  T  =  5,  NOF  =  total  number  of  frequencies  to  be  run  includ¬ 
ing  the  first  case.  After  NOF  cases  are  run,  the  program  reads 
in  the  next  complete  set  of  data. 

(4)  The  following  equations  were  used  for  P  and  Q  in  the  appro¬ 
priate  ranges 

p  +  jQ  =  - - +  j  -  .0368  + -f  n  f  — )  + 

2.54  +  1  34  r^-  ■’  [  ^  V  r^.  /  4.24J 

for  0  1  r  1  5  (101) 

-  1  ,  .  1 
2,  54  +  1  .  34  r^.  .  56  +  1  42  r^, 

for  .  5  <  r^  <  6 


(102) 
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E.  Theoretical  Antenna  Patterns 


Ground  wave  azimuth  patterns  were  calculated  and  reported  in 
Figure  3  of  Part  I  of  this  report.  These  earlier  patterns  were  displayed 
as  a  function  of  antenna  length  and  velocity  ratio,  and  were  calculated  by 
specifying  the  velocity  ratio  and  other  parameters  in  addition  to  antenna 
geometry.  Tlius  these  earlier  patterns  while  quite  valid  for  the  values 
chosen  did  not  relate  directly  to  the  earth  constants.  The  theory  has  since 
been  improved, and  the  data  are  obtained  without  assuming  velocity  ratio 
values  or  attenuation  constants.  It  is  now  only  necessary  to  specify  the 
antenna  geometry  and  the  earth  constants  (conductivity  and  dielectric 
constants).  The  computer  program  (see  Section  IV)  will  determine  all  the 
intermediate  parameters  previously  assumed  and  calculate  the  proper 
pattern,  either  azimuth  or  elevation, in  addition  to  the  effective  height. 

In  the  patterns  which  follow,  the  data  have  been  grouped  according 
to  certain  fixed  conditions  and  the  frequency  (or  some  other  characteristic) 
varied.  For  instance,  in  group  I  the  patterns  are  for  antennas  111.9  meters 
long  or  shorter  with  the  same  outer  radius.  The  earth  conductivity  has 
been  assumed  to  be  .03  mho/meter  and  the  dielectric  constant  equal  to  12 
though  any  other  reasonable  values  could  have  been  chosen.  The  chosen 
values  correspond  to  the  D/F  array  which  was  described  in  Part  I,  page  67, 
and  which  has  been  the  subject  of  extensive  experimental  investigation  as 
described  elsewhere  in  this  report. 

The  analysis  and  associated  computer  program  also  permits  cal¬ 
culation  of  theoretical  values  for  the  attenuation  constant,  velocity  ratio, 
characteristic  impedance  and  wave  tilt  angle.  These  calculated  values  are 
also  given  for  most  of  the  patterns  which  follow. 

The  computer  patterns  are  divided  into  four  groups. 

Group  I  consists  of  single  and  multiple  element  arrays,  most  111.9 
meters  long,  and  all  applicable  to  the  present  Navy  D/F  array  at  the 
Southwest  Research  Institute  D/F  tower  (outer  radius  136.65  meters). 

This  group  also  includes  some  patterns  for  antennas  and  arrays  which 
could  be  formed  by  increasing  the  inner  radius  so  as  to  shorten  the  antennas 
but  keep  the  outer  radius  fixed. 

Group  II  consists  of  similar  patterns  for  antennas  300  meters  long 
with  a  fixed  outer  array  radius  of  325  meters.  This  group  includes  some 
patterns  for  arrays  which  could  be  formed  by  increasing  the  inner  radius 
and/or  decreasing  the  outer  radius  of  the  array.  These  patterns  are 
primarily  applicable  to  the  Signal  Corps  D/F  site. 
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Croup  III  consists  of  patterns  for  pure  horizontal  polarization  and 
for  mixture  of  vertical  and  horizontal  polarization  usually  equal  amplitudes 
and  in  phase  (linear  polarization  at  45*).  These  are  for  various  lengths. 

Croup  IV  consists  of  a  number  of  miscellaneous  patterns  plotted  to 
show  certain  special  conditions,  such  as  the  elevation  angle  at  which  a 
maximum  occurs  as  a  function  of  velocity  ratio. 

The  patterns  which  appear  in  the  four  groups  are  listed  at  the  end 
of  this  section. 

Data  Croup  I  first  shows  in  Figure  I- la  through  I- 12b  both  the 
azimuth  and  elevation  patterns  from  1  to  40  mcs  for  a  single  antenna  1 
meter  above  ground.  The  "a"  figure  is  the  azimuth  pattern,  the  "b"  figure 
the  elevation  pattern.  The  source  field  is  assumed  to  be  a  vertically  polar¬ 
ized  ground  wave  for  the  azimuth  patterns,  and  a  vertically  polarized  sky 
wave  for  the  elevation  patterns.  Many  of  the  elevation  patterns  have  been 
plotted  from  0*  to  180®  elevation  angle  at  zero  azimuth,  thus  these  also 
show  the  elevation  pattern  of  the  back  lobe  at  180®  azimuth.  Both  rectangu¬ 
lar  log  plots  and  polar  plots  are  shown. 

It  is  characteristic  of  the  sky  wave  elevation  pattern  to  return  to 
zero  amplitude  at  zero  elevation.  This  is  not  a  mathematical  limitation; 
for  a  smooth  earth  it  is  quite  correct  for  the  sky  wave  component  alone 
(no  ground  wave)  and  always  results  when  the  ground  wave  or  surface  wave 
component  is  ignored.  The  relationship  between  the  sky  wave  and  ground 
wave  for  a  real  transmitter  depends  on  the  distance  to  the  transmitter, 
and  yet  there  seems  to  be  no  reason  to  pick  a  particular  transmitter  dis¬ 
tance  as  representative  of  the  manner  in  which  the  sky  wave  blends  into 
the  ground  wave.  However,  it  is  evident  that  if  a  transmitter  were  to  be 
varied  in  elevation  from  some  height  down  to  zero  degrees,  and  simul¬ 
taneously  was  adjusted  to  maintain  the  incident  field  strength  constant  at  the 
Beverage  receiving  antenna,  as  zero  elevation  was  approached,  the  eleva¬ 
tion  response  pattern  would  not  return  to  zero  at  zero  elevation.  At  zero 
elevation,  in  such  a  case,  the  response  would  be  that  typical  of  the  ground 
wave.  This  amplitude  value  (at  zero  elevation  or  the  ground  wave  case)  has 
been  indicated  on  each  polar  plot  of  the  elevation  pattern  as  a  short  line  at 
the  appropriate  amplitude,  connecting  the  zero  elevation  axis  with  the  sky 
wave  pattern. 

For  instance,  this  value  is  shown  at  the  normalized  amplitude  of  .44 
in  the  polar  plot  of  Figure  I- lb.  It  will  be  noted  that  the  transition  region 
from  sky  wave  to  ground  wave  in  all  cases  in  near  or,  at  the  higher  frequen¬ 
cies,  less  than  the  wave  tilt  angle.  Thus  the  exact  details  of  the  transition 
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at  these  low  angles  of  elevation  is  not  of  great  importance  because  no 
arriving  wave,  either  ground  wave  or  sky  wave,  ever  appears  to  arrive 
at  a  lower  elevation  angle  than  the  wave  tilt  angle. 

The  azimuth  and  elevation  patterns  show  a  decrease  in  both  beam- 
width  and  angle  of  elevation  for  maximum  response,  as  the  frequency  is 
increased.  In  general  the  back  lobe  of  the  azimuth  pattern  is  also  decreased 
with  increasing  frequency.  Azimuth  patterns  at  higher  angles  of  elevation 
are  given  later  in  data  Group  I. 

Note  that  the  elevation  patterns  show  the  main  lobe  to  be  down  to 
approximately  10*  elevation  at  30  mcs.  Since  most  long  range  sky  waves 
arrive  at  angles  below  15"  or  20"  elevation, it  is  evident  that  good  azimuth 
response  patterns  should  be  obtained  at  all  frequencies  up  to  30  mcs  for 
distant  transmitters,  say  greater  than  500  miles.  This  was  found  to  be 
the  case  during  the  bearing  accuracy  tests  reported  in  Section  III.  The 
beamwidth  and  beam  elevation  data  for  many  of  these  patterns  are  summarized 
in  Part  II  in  greater  detail 

Figures  1-13  through  1-18  show  the  change  in  the  azimuth  pattern 
as  a  function  of  the  number  of  antennas  summed  for  112-meter  antennas  2" 
apart  at  10  mcs.  It  is  evident  that  the  beamwidth  is  significantly  reduced 
even  though  no  phasing  is  used.  Figures  I- 18  through  I- 28  show  azimuth 
patterns  for  the  same  21-antenna  array,  as  the  elevation  angle  is  increased. 

It  can  be  seen  that  a  good  azimuth  pattern  is  obtained  only  for  elevations 
below  approximately  the  3-db  point  on  the  upper  side  of  the  main  beam. 

Above  this  elevation  the  azimuth  pattern  is  split.  This  result  is  the  first 
of  several  calculations  which  suggest  that  optimum  patterns  for  higher 
angles  of  elevation  will  be  obtained  with  short  elements  at  a  large  array 
radius.  Figure  1-28  is  the  elevation  pattern  for  this  example  and  shows 
that  the  poor  azimuth  patterns  obtained  around  38"  to  40"  elevation  are 
associated  with  a  minimum  in  the  elevation  pattern.  It  is  apparent  that 


26.  This  is  because  the  resultant  field,  which  is  the  sum  of  the  incident 
and  reflected  rays,  is  tilted  forward  at  an  angle  which  approaches  the  wave 
tilt  angle,  at  low  angles  of  elevation,  but  never  achieves  a  value  less  than 
the  ground  wave  tilt  angle.  For  further  reference  on  this  point,  see  Jordan, 
E.  C.,  Electromagnetic  Waves  and  Radiating  Systems,  Prentice-Hall,  New 
York,  1950,  Chapter  16;  Piggott,  W.  R.,  "A  Method  of  Determining  the 
Polar  Diagrams  of  Long  Wire  and  Horizontal  Rhombic  Aerials,  "  Department 
of  Science  and  Industrial  Research,  Radio  Research  Special  Report  No.  16, 
London,  His  Majesty's  Stationery  Office,  1948,  page  5.  See  also  Feldman, 
C.  B.,  "The  Optical  Behavior  of  the  Ground  for  Short  Radio  Waves,  "  par¬ 
ticularly  page  789,  Proceedings  of  the  IRE,  Vol.  21,  No.  6,  June  1933. 
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if  the  antenna  length  is  decreased  this  minimum  will  go  up  in  elevation 
improving  the  azimuth  pattern  at  higher  angles. 

Figures  1-29  thru  1-32  show  that  for  nonequally  spaced  arrays 
narrow  azimuth  patterns  may  be  obtained  at  the  cost  of  high  side  lobe 
levels. 

Figures  1-33  through  1-44  are  largely  a  repeat  of  Figures  I-l 
through  12  at  a  slightly  lower  height  above  ground.  The  selected  value 
of  .  85  meter  is  the  average  height  of  the  1 12  meter  long  antennas  in 
the  experimental  arrays  described  in  Part  II  of  this  report.  The  results 
are  similar  to  those  of  the  previous  figures  except  that  attenuation  is 
slightly  greater,  and  the  patterns  are  slightly  modified  in  shape. 

Figures  1-45  through  1-55  show  reduction  in  beamwidth  with  an 
increased  inner  radius  as  a  function  of  the  number  of  antennas  summed. 

Note  that  both  a  more  favorable  elevation  pattern  and  a  narrower  azimuth 
pattern  are  obtained.  Figures  1-56  through  1-66  show  the  same  result 
for  a  still  shorter  element.  Figures  1-57  through  1-82  then  explore  azi¬ 
muth  patterns  vs  elevation  angle  at  a  number  of  antennas  which  produce 
minimum  azimuth  beamwidth  at  zero  elevation. 

These  patterns  with  larger  inner  radii  and  the  same  as  previous 
outer  radius  show  much  improved  elevation  patterns  and  in  addition 
reduced  polarization  error  as  will  be  shown  in  Group  III.  Twenty-five 
meter  long  antennas  show  good  azimuth  patterns  at  all  angles  of  elevation 
of  a  summed  group  of  19  antennas  each  2®  apart  witli  an  outer  radius  of 
1 36 . 65  meters. 

The  figures  in  Group  II  are  concerned  with  antennas  300  meters 
long  or  shorter  but  on  a  site  capable  of  holding  300  meter  long  antennas. 

Figures  II- 1  through  11-12  show  single  element  patterns  vs  frequency  for 
both  azimuth  and  elevation  planes.  These  patterns  are  for  antennas  one 
meter  high.  It  is  evident  that  the  longer  element  has  produced  a  more 
favorable  azimuth  pattern  and  a  less  favorable  elevation  pattern  compared 
to  the  112  meter  long  antennas  (unless  only  low  angles  of  arrival  are  of 
interest).  Figures  II-13  thru  11-24  show  the  elevation  patterns  for  the 
same  antenna  at  a  height  of  two  meters. 

Figures  11-25  thru  11-32  show  azimuth  patterns  vs  elevation  at 
6.0  mes  for  a  single  element  2  meters  high  and  300  meters  long.  Again, 
as  with  the  112-meter  antennas,  split  azimuth  patterns  are  obtained  at 
elevations  above  the  main  lobe.  Figures  11-33  thru  11-36  show  variation 
with  frequency  between  2.0  and  6.0  mes  at  an  elevation  of  20®. 

Figures  11-37  thru  11-42  show  beamwidth  reduction  with  increased 
number  of  summed  antennas.  Figures  11-43  thru  11-52  show  an  increased 
beamwidth  reduction  for  a  greater  inner  radius  and  less  outer  radius. 
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Figures  11-53  thru  11-64  show  patterns  for  the  same  length  previously 
shown  with  an  outer  radius  of  325  meters.  Minimum  beamwidth  of 
approximately  12*  is  obtained  with  15  antennas  at  10  mcs. 

Figures  11-65  to  11-70  show  azimuth  patterns  vs  number  of  antennas 
for  a  shorter  element  length  than  the  previous  set,  but  the  same  outer 
radius.  Figures  11-71  thru  74  show  a  single  element  pattern  variation  with 
height  above  ground. 

Figures  11-75  thru  11-78  show  split  sector  arrays  of  antennas  all 
2®  apart  and  illustrate  how  still  narrower  azimuth  patterns  may  be  obtained 
at  the  cost  of  higher  side  lobe  level.  A  pattern  designated  7  on  1  off  7  on 
means  14  antennas  in  two  sectors  of  seven  each.  The  two  sectors  are 
separated  by  a  gap  of  4  degrees  (one  missing  antenna).  Similarly  the 
6  on  3  off  6  on  array  means  two  6  element  sectors  separated  by  an  8“ 
gap  (three  missing  antennas).  These  type  arrays  can  be  obtained  with 
the  present  commutator. 

Further  details  and  discussions  of  the  various  data  obtained  from 
the  patterns  are  given  in  Part  II  of  this  report. 

Patterns  in  Group  III  are  concerned  with  the  effects  of  horizontally 
polarized  fields,  or  elliptically  polarized  fields.  The  first  20  patterns 
shown  in  Figures  III-l  thru  III-20  show  the  response  of  one  and  two 
meter  high  single  elements  to  pure  horizontal  polarization.  These  patterns 
are  symmetrical  and  have  a  null  in  the  direction  of  the  antenna  line. 

Patterns  like  those  in  the  first  20  figures  of  this  group  have  rarely 
been  observed  experimentally  because  the  horizontally  polarized  compo¬ 
nents  almost  always  occur  in  the  presence  of  vertically  polarized  com¬ 
ponents.  In  general  the  response  of  the  antenna  is  greater  to  vertically 
polarized  than  horizontally  polarized  fields  of  equal  amplitude.  This  is 
because  only  the  vertically  polarized  field  occurs  in  a  manner  favorable 
to  a  traveling  wave  buildup  on  the  antenna.  Figures  III-21  thru  III-24 
illustrate  this  for  low  angles  of  elevation  on  longer  elements. 

Figures  III-21  thru  III-24  show  that  polarization  error  is  small  at 
low  angles  of  elevation  but  may  increase  above  the  main  beam.  It  can 
also  be  seen  that  the  patterns  are  increasingly  asymmetrical  as  the  elevation 
angle  increases  (or  the  effect  of  the  horizontally  polarized  component 
increases).  This  asymmetry  also  appears  in  practice  when  signals  are 
observed  which  are  changing  in  polarization.  From  the  calculations  it 
would  be  expected  that  a  sky  wave  signal  would  appear  to  fade  and  slightly 
move  in  azimuth  for  low  angles  and  long  antennas.  If  the  polarization 
error  is  slight,  the  pattern  would  appear  to  fade  on  the  bearing  indicated 
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when  good  symmetry  occurred.  This  effect  has  been  observed  in 
practice  and  generally  the  fading  drops  the  amplitude  of  the  displayed 
pattern  well  below  the  maximum  before  a  noticable  polarization  shift 
occurs.  One  would  expect  then  on  occasion  to  observe  pure  horizontally 
polarized  patterns  at  the  lowest  part  of  some  fades;  that  is,  patterns  such 
as  those  in  Figures  III-l  thru  III-20.  This  effect  has  also  been  observed 
as  expected  although  not  frequently.  The  fact  that  these  expected  patterns 
have  been  observed  however,  lends  support  to  the  general  method  of 
analysis  given  in  this  report. 

If  the  analysis  given  is  correct,  then  for  long  antennas  at  low 
angles,  polarization  error  will  not  be  appreciable  for  two  reasons'  (1)  its 
measured  value  will  be  low,  and  (2)  the  appearance  of  the  pattern  will  be 
such  that  it  can  be  recognized  that  at  a  given  instant  some  polarization 
error  exists,  i.e.,  asymmetry  and  slight  back  and  forth  swinging  as  the 
polarization  shifts,  so  that  bearings  can  be  read  at  an  optimum  time. 

In  Part  III  of  this  report  results  are  presented  for  operators  with  no 
experience  at  pattern  interpretation  vs  those  using  symmetry  as  criteria. 
Results  bear  out  what  has  been  stated  above. 

Further  consideration  suggests  that  at  high  angles  of  elevation  the 
polarization  error  should  increase  if  the  antenna  element  is  long  (Figure 
III-23).  This  leads  to  the  same  design  direction  indicated  earlier,  that 
of  using  shorter  elements  to  improve  the  elevation  pattern. 

This  is  explored  in  Fi.'nrcs  III-30  thru  III-33  where  the  azimuth 
patterns  for  an  array  of  19  antennas  each  25  meters  long  are  shown  as  a 
function  of  angle  of  elevation.  It  is  evident  that  an  improvement  has  been 
obtained  in  that  the  polarization  error  is  small  up  to  30°.  The  previous 
remarks  concerning  asymmetry  or  a  criterion  for  error  detection  are 
equally  applicable  in  this  case.  Further  and  more  detailed  consideration 
of  response  to  horizontal  polarization  is  given  in  Part  II  of  this  report. 

Patterns  shown  by  Figures  IV- 1  thru  IV-5  are  a  miscellaneous 
group  included  to  illustrate  two  methods  of  pattern  improvement  which 
have  not  yet  been  attempted  experimentally.  Figures  IV- 1  thru  IV-4  show 
the  variation  in  a  typical  long  element  elevation  pattern  with  artificial 
variation  of  the  velocity  ratio  on  the  antenna  As  the  velocity  ratio  is 
increased  to  values  greater  than  one,  the  angle  of  elevation  of  the  main 
beam  is  increased.  Thus  this  could  be  a  method  for  controlling  response 
for  short  range  signal  arriving  at  high  angles.  (A  velocity  ratio  greater 
than  one  implies  a  phase  velocity  faster  than  the  velocity  of  light  This 
may  be  achieved  for  instance  by  inserting  series  capacitance  in  the  antenna 
as  was  done  by  Beverage.  9) 


Figure  IV-5  shows  a  method  whereby  a  very  sharp  azimuth  beam 
may  be  obtained  without  using  a  phased  array.  The  patterns  were  obtained 
by  reversing  the  phase  of  a  pair  of  antennas  in  the  array  before  summing. 
The  null  between  the  main  beam  and  the  first  side  lobe  occurs  in  the 
vicinity  of  the  direction  where  the  phase  reversal  occurs  in  azimuth.  An 
array  of  180  antennas  could  be  used  to  produce  such  patterns  without  a 
phasing  goniometer  by  using  multicouplers  (or  other  devices  with  both  0* 
and  180*  phase  outputs)  and  a  360-input  commutator  (2  inputs  per  antenna). 


GROUP  I  -  112-METER  ANTENNAS 


Figure 

Pattern 

Freq. 

No. 

Plane 

in  Mcs 

I-la 

Azimuth 

l.O 

I-lb 

Elevation 

l.O 

I-2a 

Azimuth 

1.5 

I-2b 

Elevation 

1.5 

I-3a 

Azimuth 

2.0 

I-3b 

Elevation 

2.0 

I-4a 

Azimuth 

3.0 

I-4b 

Elevation 

3.0 

I-5a 

Azimuth 

4.0 

All  single  antennas  one  meter  high. 

I-5b 

Elevation 

4.  0 

Azimuth  patterns  ground  wave 

I-6a 

Azimuth 

6.0 

>.  vertical  polarization,  elevation 

I-6b 

Elevation 

6.0  1 

patterns  sky 

wave  vertical  polari- 

I-7a 

Azimuth 

8.0 

zation  at  zero  azimuth. 

I-7b 

Elevation 

8.0 

I-8a 

Azimuth 

10.0 

I-8b 

Elevation 

10.0 

I-9a 

Azimuth 

15.0 

I-9b 

Elevation 

15.0 

I-lOa 

Azimuth 

20. 0 

I-lOb 

Elevation 

20.0 

I-lla 

Azimuth 

30.0 

I-llb 

Elevation 

30.0 

I-12a 

Azimuth 

40.0 

I-12b 

Elevation 

40.0 

1-13 

Azimuth 

10.0 

Single  Antenna 

1-14 

Azimuth 

10.0 

3  Antennas 

All  antennas  2°  apart, 

1-15 

Azimuth 

10.0 

5  Antennas 

1  meter  high,  zero 

1-16 

Azimuth 

10.0 

9  Antennas  j 

f  elevation 

1-17 

Azimuth 

10.0 

1 5  Antennas  J 

1-18 

Azimuth 

10.0 

Zero  Elevationi 

1-19 

Azimuth 

10.0 

10°  Elevation 

1-20 

Azimuth 

10.0 

20°  Elevation 

1-21 

Azimuth 

10.0 

30°  Elevation 

All  patterns 

1-22 

Azimuth 

10.0 

32°  Elevation 

21  antennas  all 

1-23 

Azimuth 

10.0 

34°  Elevation 

>■  2°  apart,  1  meter 

1-24 

Azimuth 

10.0 

36°  Elevation 

high,  vertical 

1-25 

Azimuth 

10.0 

38°  E\evation 

polarization 

1-26 

Azimuth 

10.0 

40°  Elevation 

1-27 

Azimuth 

10.0 

42°  Elevation 

1-28 

Elevation 

10.0 

Zero  Azimuth 
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CROUP  I  -  i  12-METER  ANTENNAS  (Cont'd) 


Figure 

Pattern 

Freq. 

No. 

Plane 

in  Mcs 

1-29 

Azimuth 

10.0 

2  Antennas,  ±20 

*  azimuth 

1-30 

Azimuth 

10.0 

6  Antennas,  ±18 

*,  ±20*,  ±22"  azimuth 

1-31 

Azimuth 

10.0 

6  Antennas,  ±16 

*,  ±20*,  ±24*  azimuth 

1-32 

Azimuth 

10.0 

6  Antennas,  ±16 

*,  ±18*,  ±20*  azimuth 

1-33 

Azimuth 

1.0 

1-34 

Azimuth 

1.5 

1-35 

Azimuth 

2.0 

1-36 

Azimuth 

3.0 

1-37 

Azimuth 

4.0 

1-38 

Azimuth 

6.0 

All  single  antennas  .  85  meter  high. 

1-39 

Azimuth 

8.0 

r  Azimuth  pattern 

s  have  vertically 

1-40 

Azimuth 

10.0 

polarized  ground  wave. 

1-41 

Azimuth 

15.0 

1-42 

Azimuth 

20.0 

1-43 

Azimuth 

30.0 

1-44 

Azimuth 

40.0 

J 

1-45 

Azimuth 

Single  Antenna 

All  antennas  2* 

1-45 

Azimuth 

3  Antennas 

apart,  one  meter 

1-47 

Azimuth 

5  Antennas 

high,  zero  eleva¬ 

1-48 

Azimuth 

7  Antennas 

tion.  Antennas  50 

1-49 

Azimuth 

9  Antennas 

meters  long,  inner 

1-50 

Azimuth 

1 1  Antennas 

radius  86  65 

1-51 

Azimuth 

1  5  Antennas 

meters,  outer 

1-52 

Azimuth 

19  Antennas 

radius  136  65 

1-53 

Azimuth 

21  Antennas 

meters.  All 

1-54 

Elevation 

Single  Antenna 

patterns  at 

1-55 

Elevation 

21  Antennas 

10  me. 

1-56 

Azimuth 

Single  Antenna 

All  antennas  2* 

1-57 

Azimuth 

3  Antennas 

apart,  one  meter 

1-58 

Azimuth 

5  Antennas 

high,  zero  eleva¬ 

1-59 

Azimuth 

7  Antennas 

tion.  Antennas  25 

1-60 

Azimuth 

9  Antennas 

meters  long,  inner 

1-61 

Azimuth 

1 1  Antennas 

>  radius  111.65 

1-62 

Azimuth 

1  5  Antennas 

meters,  outer 

1-63 

Azimuth 

1  7  Antennas 

radius  136. 65 

1-64 

Azimuth 

21  Antennas 

meters.  All 

1-65 

Elevation 

Single  Antenna 

patterns  at 

1-66 

Elevation 

21  Antennas 

J 

10  me. 

65 


CROUP  I  -  112-METER  ANTENNAS  (Cont’d) 


Figure 

Pattern 

Freq. 

No. 

Plane 

in  Mcs 

1-67 

Azimuth 

5*  Elevation 

All  antennas  2* 

1-68 

Azimuth 

10*  Elevation 

apart,  one  meter 

1-69 

Azimuth 

20*  Elevation 

high.  Antennas 

1-70 

Azimuth 

25*  Elevation 

25  meters  long, 

1-71 

Azimuth 

30*  Elevation 

inner  radius 

1-72 

Azimuth 

35*  Elevation 

111.6  meters. 

1-73 

Azimuth 

40*  Elevation 

outer  radius 

1-74 

Azimuth 

45*  Elevation 

^  1 36.  6  meters. 

1-75 

Azimuth 

50*  Elevation 

10  me.  Verti¬ 

1-76 

Azimuth 

55*  Elevation 

cally  polarized 

1-77 

Azimuth 

60*  Elevation 

sky  wave.  Azi¬ 

1-78 

Azimuth 

65*  Elevation 

muth  plane 

1-79 

Azimuth 

70*  Elevation 

patterns  for 

1-80 

Azimuth 

75*  Elevation 

sum  of  19 

1-81 

Azimuth 

80*  Elevation 

antennas. 

1-82 

Azimuth 

85*  Elevation 
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GROUP  II  -  300-METER  ANTENNAS 


Figure 

Pattern 

Freq. 

No. 

Plane 

in  Mcs 

Il-la 

Azimuth 

1.0 

Il-lb 

Elevation 

1.0 

II-2a 

Azimuth 

1.5 

II-2b 

Elevation 

1.5 

II-3a 

Azimuth 

2.0 

II-3b 

Elevation 

2.0 

II -4a 

Azimuth 

3.0 

II-4b 

Elevation 

3.0 

II-5a 

Azimuth 

4.0 

All  single  antennas  one  meter 

II-5b 

Elevation 

4.0 

high.  Azimuth  patterns  ground 

II -6a 

Azimuth 

6.0 

wave  vertical  polarization,  ele¬ 

II-6b 

Elevation 

6.0 

r 

vation  patterns  sky  wave  vertical 

II-7a 

Azimuth 

8.0 

polarization  at  zero  azimuth. 

II -7b 

Elevation 

8.0 

II-8a 

Azimuth 

10.0 

II -8b 

Elevation 

10.0 

II-9a 

Azimuth 

15.0 

II-9b 

Elevation 

15.0 

II- 10a 

Azimuth 

20.0 

Il-lOb 

Elevation 

20.0 

Il-lla 

Azimuth 

30.0 

Il-llb 

Elevation 

30.0 

II-12a 

Azimuth 

40.0 

II-12b 

Elevation 

40.0  ^ 

11-13 

Elevation 

1.0  ^ 

11-14 

Elevation 

1.5 

11-15 

Elevation 

2.0 

11-16 

Elevation 

3.0 

11-17 

Elevation 

4.0 

All  single  antennas  two  meters 

11-18 

Elevation 

6.0 

> 

high.  Sky  wave,  vertical 

11-19 

Elevation 

8.0 

polarization,  at  zero  azimuth. 

11-20 

Elevation 

10.0 

11-21 

Elevation 

15.0 

11-22 

Elevation 

20. 0 

11-23 

Elevation 

30.0 

11-24 

Elevation 

40.0 
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GROUP  II  -  300-METER  ANTENNAS  (Cont'd) 


Figure 

Pattern 

Freq. 

No. 

Plane 

in  Mcs. 

11-25 

Azimuth 

6.0 

Zero  Elevation 

'I  All  single 

11-26 

Azimuth 

6.0 

10*  Elevation 

antennas,  two 

11-27 

Azimuth 

6.0 

20*  Elevation 

meters  high. 

11-28 

Azimuth 

6.0 

30*  Elevation 

sky  wave  verti¬ 

11-29 

Azimuth 

6.0 

40*  Elevation 

cal  polarization. 

11-30 

Azimuth 

6.0 

50*  Elevation 

Antennas  300 

11-31 

Azimuth 

6.0 

60*  Elevation 

meters  long. 

11-32 

Elevation 

6.0 

Zero  Azimuth 

> 

11-33 

Azimuth 

2.0  1 

All  single  antennas,  two  meters 

11-34 

Azimuth 

3.0  1 

high,  elevation  angle  20*,  sky 

11-35 

Azimuth 

4.0  ( 

wave  vertical  polarization. 

11-36 

Azimuth 

6.0  J 

Antenna  length  300  meters. 

11-37 

Azimuth 

Single  Antenna' 

All  antennas  2* 

11-38 

Azimuth 

3  Antennas 

apart,  two  meters 

11-39 

Azimuth 

5  Antennas 

high,  sky  wave 

11-40 

Azimuth 

9  Antennas 

vertical  polariza¬ 

II-41 

Azimuth 

15  Antennas 

tion.  Antennas  300 

11-42 

Azimuth 

21  Antennas 

meters  long,  inner 

radius  25  meters, 
outer  radius  325 
meters,  frequency 

10  me. 

11-43 

Azimuth 

Single  Antenna' 

All  antennas  2® 

11-44 

Azimuth 

3  Antennas 

apart,  one  meter 

11-45 

Azimuth 

5  Antennas 

high,  sky  wave 

11-46 

Azimuth 

7  Antennas 

vertical  polariza- 

11-47 

Azimuth 

9  Antennas 

.  tion.  Antennas  100 

11-48 

Azimuth 

1 1  Antennas 

meters  long,  inner 

11-49 

Azimuth 

13  Antennas 

radius  125  meters, 

11-50 

Azimuth 

15  Antennas 

outer  radius  225 

11-51 

Azimuth 

21  Antennas 

meters,  frequency 

11-52 

Elevation 

Single  Antenna 

10  me. 

68 


CROUP  11  -  300-METER  ANTENNAS  (Cont'd) 


Figure 

Pattern 

Freq. 

• 

No. 

Plane 

in  Mcs 

11-53 

Azimuth 

Single  Antenna 

- 

11-54 

Azimuth 

3  Antennas 

11-55 

Azimuth 

5  Antennas 

All  antennas  2*  apart, 

11-46 

Azimuth 

7  Antennas 

one  meter  high,  ver¬ 

11-57 

Azimuth 

9  Antennas 

tical  polarization. 

11-58 

Azimuth 

1 1  Antennas 

^  Antennas  100  meters 

11-59 

Azimuth 

13  Antennas 

long,  inner  radius 

11-60 

Azimuth 

15  Antennas 

225  meters,  outer 

11-61 

Azimuth 

19  Antennas 

radius  325  meters 

11-62 

Azimuth 

21  Antennas 

frequency  10  mcs. 

11-63 

Elevation 

Single  Antenna 

11-64 

Elevation 

15  Antennas 

11-65 

Azimuth 

3  Antennas 

All  antennas  2®  apart. 

11-66 

Azimuth 

5  Antennas 

one  meter  high,  ver- 

11-67 

Azimuth 

7  Antennas 

^  tical  polarization. 

11-68 

Azimuth 

9  Antennas 

Antennas  50  meters 

11-69 

Azimuth 

1 1  Antennas 

long,  inner  radius 

11-70 

Azimuth 

15  Antennas 

275  meters,  outer 

radius  325  meters. 

frequency  10  mcs. 

11-71 

Azimuth 

15  Meters  high 

Single  antenna,  vertical 

11-72 

Azimuth 

1  Meter  high 

polarization,  antenna 

11-73 

Azimuth 

2  Meters  high 

300  meters  long,  fre¬ 

11-74 

Azimuth 

4  Meters  high 

quency  10  mcs,  eleva¬ 

> 

tion  angle  20® . 

11-75 

Azimuth 

"7  on,  1  off,  7  on" 

15  antennas  in  a  split 

11-76 

Azimuth 

"6  on,  3  off,  6  on" 

sector  array,  all  2® 

11-77 

Azimuth 

"5  on,  5  off,  5  on" 

apart  whether  on  or 

11-78 

Azimuth 

"4  on,  7  off,  4  on" 

off.  Seven  on,  1  off. 

I  7  on  means  7  antennas 
^  2®  apart  (a  12®  sector) 
then  a  4®  sector  with 
no  antennas,  then  7 
antennas  2®  apart  again 
to  complete  symmetry. 
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CROUP  III  -  HORIZONTAL  POLARIZATION 


Figure 

Pattern 

Freq. 

No. 

Plane 

m  Mcs 

III-l 

Azimuth 

1.0 

max  " 

All  single  antennas,  pure  horizontal 

III -2 

Azimuth 

1.5 

max 

polarization,  one  meter  high.  Maxi¬ 

III-3 

Azimuth 

2.0 

max 

mum  and  minimum  refers  to  eleva¬ 

III-4 

Azimuth 

2.0 

min 

tion  angle  chosen.  Maximum  is 

III-5 

Azimuth 

3.0 

max 

elevation  angle  of  a  maximum  on 

III-6 

Azimuth 

3.0 

min 

elevation  plane  for  sky  wave  vertical 

III-7 

Azimuth 

4.0 

max 

polarization.  Minimum  is  elevation 

III-8 

Azimuth 

4.0 

min 

angle  of  a  minimum  on  elevation 

III-9 

Azimuth 

6.0 

max 

plane  pattern  for  sky  wave  vertical 

III-IO 

Azimuth 

6.0 

min 

polarization. 

III-ll 

Azimuth 

1.0 

max 

All  single  antennas,  pure  horizontal 

III-12 

Azimuth 

1.5 

min 

polarization,  two  meters  high.  Maxi¬ 

III-13 

Azimuth 

2.0 

max 

mum  and  minimum  refers  to  eleva¬ 

III-14 

Azimuth 

2.0 

min 

tion  angle  chosen.  Maximum  is 

III-15 

Azimuth 

3.0 

max 

y  elevation  angle  of  a  maximum  on 

III-16 

Azimuth 

3.0 

min 

elevation  plane  for  sky  wave  vertical 

III-17 

Azimuth 

4.0 

max 

polarization.  Minimum  is  elevation 

III- 18 

Azimuth 

4.0 

min 

angle  of  a  minimum  on  elevation 

III-19 

Azimuth 

6.0 

max 

plane  pattern  for  sky  wave  vertical 

III-20 

Azimuth 

6.0 

min 

J 

polarization. 

III-21 

Azimuth 

6.0 

Single 

Antenna,  10°  Elevation  Mixed  polari- 

III-22 

Azimuth 

6.0 

Single 

Antenna,  20°  Elevation  zation,  vertical 

III-23 

Azimuth 

6.0 

Single 

Antenna,  30°  Elevation  plus  horizon- 

III-24 

Azimuth 

6.0 

Single 

Antenna,  40°  Elevation  tal 

III-25 

Azimuth 

6.0 

.5  meter  high  Single  antennaS; 

III-26 

Azimuth 

6.0 

1.0  meter  high  20°  elevation, 

III-27 

Azimuth 

6.0 

2.0  meters  high  shows  both 

III-28 

Azimuth 

6.0 

4.0  meters  high  vertical  and 

III-29 

Azimuth 

10.0 

shows 

polarization  error  vs  mixed  polari- 

number  of  antennas  zation  patterns. 

III-30 

Azimuth 

All  azimuth  plane  patterns  at  10  me. 

III-31 

Azimuth 

Sum  of  19  antennas,  2°  apart,  one 

III-32 

Azimuth 

meter  high.  Antennas  25  meters  long, 

III-33 

Azimuth 

inner  radius  111.6  meters,  outer 

III-34 

Azimuth 

radius  136.6  meters.  Linear  polariza- 

tion  at  45°  in  phase  (equal  in  phase  com¬ 
ponents  of  vertical  and  horizontal 
polarization).  10  mcs. 
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GROUP  IV 


Figure 

Pattern 

Freq. 

No. 

Plane 

in  Mcs 

IV-1 

Elevation 

10.0 

Velocity  Ratio  =1.00 

Single  antennas 

IV-2 

Elevation 

10.0 

Velocity  Ratio  =  1.  10 

111.9  meters 

IV-3 

Elevation 

10.0 

Velocity  Ratio  =  1.20 

long,  one  meter 

IV-4 

Elevation 

10.0 

*  Velocity  Ratio  =  1.  30 

high. 

IV-5 

Azimuth 

10.0 

Various  patterns  for  group  of  four 

antennas,  2  pairs  180“ 

out  of  phase. 

j 


mm 


' 
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i 
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F.  Comments  and  Conclusions 


Theoretical  characteristics  of  a  circular  array  of  radially 
arranged  Beverage  antennas  have  been  investigate  1  with  respect  to 
direction  finding  applications  in  the  medium,  high  frequency  and  very 
high  frequency  ranges.  Results  suggest  that  the  Beverage  antenna  may 
be  used  as  an  effective  antenna  element  in  wide  aperture  circular  direc¬ 
tion  finding  arrays  over  wide  frequency  ranges  normally  associated 
with  frequency  independent  antennas 

The  Beverage  antenna  consists  simply  of  a  long  horizontal  wire, 
parallel  to  the  surface  of  the  earth  and  at  a  low  height  above  ground. 

It  is  always  arranged  radially  from  the  receiving  point  to  the  direction 
of  desired  reception  since  it  is  an  end  fire  antenna  operating  in  a  traveling 
wave  mode  .  A  typical  element  for  the  2  to  30  me  range  may  be  con¬ 
structed  for  less  than  $50  (the  cost  of materials)and  it  is  therefore 
possible  to  construct  a  circular  array  of  180  such  antennas  on  a  900-ft 
diameter  site  for  less  than  a  total  cost  of  $20,000  exclusive  of  receiving  and 
commutating  equipment,  and  input  cables.  In  view  of  such  an  unusually 
low  cost  for  a  wide  aperture  antenna  system,  covering  more  than  the 
entire  high  frequency  range  ,  it  is  of  interest  to  make  a  detailed  theo¬ 
retical  and  experimental  examination  of  performance  characteristics 
Both  investigations  have  been  under  way  at  this  laboratory  for  about 
four  years  A  new  and  detailed  theoretical  investigation  has  been  under 
way  for  the  past  two  years,  and  a  number  of  broad  conclusions  may  now 
be  stated  This  report  (Part  IV)  is  almost  entirely  concerned  with  theory 
except  where  experimental  work  has  provided  confirmation  or  where 
directions  of  theoretical  investigation  have  been  indicated  Detailed 
conclusions  are  given  in  Part  II  and  Part  III  of  this  report  which  concern 
theoretical  results  as  compared  with  measured  results  on  existing  system.s 

1  Accuracy  Expected  Based  on  Theory 

It  was  reported  in  1962  that  low  standard  deviations  could 
be  obtained  on  long  range  randomly  polarized  sky  wave  signals  in  the 
frequency  range  of  3  to  30  me  In  1961  a  figure  of  1  “  was  measured  for 
36  antennas  300  meters  long  looking  in  a  72®  sector  toward  the  northeast 
quadrant  In  1961  a  standard  deviation  of  4  28°  was  also  obtained  with  36 
112-meter-long  antennas  covering  360°  during  the  same  testing  period 
Further  accuracy  experiments  have  been  conducted  since  that  time  showing 
improved  results  (see  Part  III)  and  some  accuracy  calculations  reported 
herein  have  been  made  for  assumed  incident  fields  of  mixed  polarization. 
Polarization  error  calculations  reported  herein  indicate  that  long  range 
horizontally  polarized  sky  wave  components  may  momentarily  introduce 
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small  errors,  up  to  a  few  degrees,  under  the  worst  conditions.  The 
appearance  of  the  patterns  when  horizontally  polarized  components  are 
present  is  slightly  asymmetrical  (or  fully  symmetrical  when  only 
vertical  polarization  is  present).  Thus  even  though  the  polarization 
error  is  low,  a  clear  criterion  exists  for  observing  when  polarization 
error  is  at  its  lowest.  The  polarization  error  does  not  appear  to  be 
significant,  and  the  results  of  experiments  reported  in  Section  III 
bear  this  out. 

2 .  Bandwidth 


The  antenna  is  essentially  a  matched  transmission  line 
and  has  an  impedance  which  is  approximately  independent  of  frequency 
Calculations  indicate  that  very  favorable  impedances  can  be  obtained 
for  a  40  to  1  frequency  range,  say,  from  1  to  40  me.  There  is  no  theo¬ 
retical  reason  why  a  constant  impedance  may  not  be  obtained  for  a  very 
wide  frequency  range,  say,  1  to  100  me  or  a  range  which  is  so  wide  as 
to  be  limited  by  other  characteristics. 

It  is  possible  to  obtain  narrow  beamwidths  in  the  azimuth 
plane  over  a  wide  frequency  range  without  phasing  the  antennas  Further¬ 
more,  it  is  possible  to  vary  the  beamwidth  quickly  at  a  given  frequency 
It  is  also  possible  to  obtain  very  favorable  elevation  patterns  so  that 
narrow  beam  azimuth  patterns  are  obtained  up  to  angles  of  elevation  of 
70®  over  the  entire  HF  frequency  range.  It  is  to  be  emphasized  that  these 
patterns  are  obtained  without  phasing  the  elements  but  are  obtained  by 
simple  summation  of  a  sector  of  antennas  as  they  occur  in  the  circular 
array  (by  commutating  the  summed  sectors,  the  patterns  are  displayed 
on  a  CRT)  Further  details  and  exact  beamwidths  which  can  be  obtained 
are  reported  in  Part  II. 

3 .  Effective  Height  and  Sensitivity 

Analysis  shows  that  the  effective  height  of  an  ideal  Beverage 
antenna  is  approximately  6f  /2  where  6  is  the  wave  tilt  in  radians  typically 
on  the  order  of  10  or  15  degrees,  and  i  is  the  antenna  length  in  meters. 
This  means  that  a  single  element  100  meters  long  may  have  an  effective 
height  on  the  order  of  five  meters  at  a  typical  frequency.  Summation  of 
arrays  of  10  to  20  elements  may  produce  effective  heights  5  to  10  times 
the  element  effective  height.  The  antennas  are,  therefore,  comparable 
in  sensitivity  to  other  wide  aperture  typical  antenna  arrays  used  for  high 
frequency  direction  finding. 
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4 .  Computer  Analysis 

A  digital  computer  program  has  been  prepared  to  make 
exact  calculations  of  azimuth  and  elevation  patterns  for  circular  arrays 
of  Beverage  antennas.  Only  a  few  hundred  patterns  have  been  computed 
at  this  time,  and  the  detailed  investigation  of  all  the  parameters  which 
may  be  varied  to  produce  optimum  results  has  only  been  initiated. 
Nevertheless  favorable  patterns  in  both  azimuth  and  elevation  planes 
for  both  vertical  and  horizontal  polarization  have  been  obtained  Design 
suggestions  for  very  narrow  beam  azimuth  patterns  with  very  broad 
elevation  patterns  are  simple  and  promising.  The  most  significant  of 
these  is  to  use  shorter  elements  with  a  larger  inner  radius  for  a  fixed 
overall  array  diameter. 

5 .  Continued  Investigation 

On  the  basis  oi  the  theoretical  investigations  conducted  so 
far,  it  is  clearly  apparent  that  the  theoretical  investigation  should  be 
continued  to  optimize  azimuth  and  elevation  patterns,  and  an  extensive 
experimental  program  should  be  initiated  to  compare  the  performance 
of  a  low  cost  Beverage  system  with  other  wide  aperture  high  frequency 
direction  finders.  It  is  suggested  that  the  present  existing  antenna 
arrays  at  Southwest  Research  Institute  be  used  for  such  tests 


G. 


Recommendations 


It  is  recommended  that  investigation  of  the  circular  Beverage 
array  for  wide  aperture  high  frequency  direction  finding  be  continued, 
utilizing  the  antennas  in  their  simple  low  cost  form.  Comparsion  with 
other  wide  aperture  high  frequency  direction  finders  should  be  made  on 
a  quantitative  basis.  Several  types  of  high  frequency  direction  finders 
are  in  use  or  have  been  suggested  and  should  be  compared.  These 
include  the  Wullenweber,  the  Doppler  array,  the  Beverage  array, 
the  interferometer  array  and  others,  and  finally  the  recently  suggested 
variation  on  the  Luneberg  Lens  Each  of  these  antenna  systems  may 
have  an  area  of  application  where  its  performance  is  near  a  maximum 
Since  the  Beverage  system  is  unquestionably  the  lowest  in  cost  of  these 
four,  and  there  is  every  indication  that  its  accuracy  performance  is 
comparable,  its  performance  should  be  measured  in  considerable  detail 


=!‘Tanner,  R.  L  ,  et.  al.  ,  "An  Investigation  of  the  HF  Wire-Grid 
Lens  Antenna  for  Direction  Finding , "  RADC- TDR-63-21 1 ,  December , 
1963,  Contract  No.  AF30(602)-2742,  TRG-West,  A  Division  of  TRG, 
Incorporated,  Menlo  Park,  California. 
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The  very  large  reduction  in  mutual  impedance,  due  to  the  current 
in  the  rail,  is  striking. 
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REPRINT  OF  "PROPAGATION  OF  HIGH-FREQUENCY  CURRENTS 
IN  GROUND  RETURN  CIRCUITS"  BY  W.  H,  WISE 


Bell  Telephone  System 


A  Derivation  of  the 

iDANCE  FORMOIA  WITHOUT  THE  ASSUMPTION 
■REQUENCY  IS  SO  LOW  THAT  POLARIZATION 
CURRENTS  MAY  BE  NEGLECTED 
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APPENDIX  III 


LETTER  FROM  VV  H  WISE  CORRECTING 
PREVIOUS  PAPERS 


94  Paterson  Road 
Fanwood,  N.  J. 
October  14,  1963 


Mr.  D.  N,  Travers 
8500  Culebra  Road 
San  Antonio  6,  Texas 

Dear  Mr.  Travers* 

I  have  mailed  some  reprints  to  you  None  of  these  monographs  con¬ 
tains  an  alternative  to  calculating  Carson's  eight  series,  modified  for  high 
frequencies. 

If  you  have  been  looking  at  Carson's  BSTJ  paper,  or  Monograph 
B-219.  you  may  have  noticed  that  his  curve  labeled  tt/Z  is  out  of  place  in  the 
neighborhood  of  r  =  2.5  in  Figure  2.  The  numerical  information  for  this 
figure  was  correctly  computed  and  was  then  given  to  a  draughtsman  who 
prepared  the  figure.  Seeing  that  the  other  curves  seemed  to  pass  through 
a  single  point  at  r  =  2  he  leaped  to  the  conclusion  that  the  tt/2  curve  should 
pass  through  this  point.  So  he  passed  the  tt/?.  curve  through  this  point, 
without  mentioning  what  he  had  done.  The  crime  was  discovered  after  the 
paper  was  out.  Later,  when  reprints  had  become  a  scarce  item,  IT&T 
asked  for  a  copy  from  which  they  could  make  copies  for  their  engineers. 

I  was  given  the  task  of  correcting  Figure  2  by  scratching  out  the  displaced 
part  of  the  curve  and  then  inking  it  in  where  it  belongs  , 

There  are  two  misprints  in. Monograph  B-1558  There  are  no  mis¬ 
prints  in  the  BSTJ  article  The  contractor  for  the  reprints  used  a  proof 
submitted  for  an  estimate  of  cost  to  make  the  reprints,  instead  of  waiting 
for  the  corrected  page  proof. 

I  find  in  my  copy  of  Monograph  B-791  a  reference  to  a  paper  by 
Sidney  Weinbaum  in  the  Journal  of  Applied  Physics,  December  1944,  Vol 

e  X"  + 

15,  No.  12,  page  840,  He  shows  how  to  evaluate  I  - =—  dx 

0  V  x2  +  p2 

by  deriving  a  differential  equation  for  it  and  then  getting  a  series  solution 
for  the  differential  equation,  The  variable  is  -jkp.  This  series  is  said  to 
be  an  improvement  over  the  older  one. 

I  do  not  know  that  anyone  has  attempted  to  improve  on  Carson's 

series . 


Yours  truly, 


/s/  W  Howard  Wise 
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APPENDIX  IV 

REPRINT  CF  "POTENTIAL  COEFFICIENTS  FOR  GROUND 
RETURN  CIRCUITS"  BY  W  H  WISE 
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